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1 Introduction 

Nanotechnology and Nanoscience are inducing a turning point in fields from 
Condensed Matter Physics and Materials Science to Chemistry and Biology. 
This is due to the new types of behaviors and properties displayed by nanos- 
tructures, ranging across traditional disciplines [83]. Interest in these systems 
is also triggered by the possibility of characterizing them morphologically 
by scanning probe microscopy techniques (SPM) [11] and by the advance in 
surface structure determination methods. Within this more specific context, 
there is an increasing need for the development of techniques and methods 
for the patterning of materials at the nanoscale. Moreover, this goal must be 
achieved in an efficient, fast and low-cost manner for eventual technological 
applications to be compatible with mass-production. In addition, the pro- 
cesses involved should allow control over the size, shape and composition of 
the nanostructures produced. 

Within this trend, there is also a great interest in developing methods for 
producing confined nanostructures, which would display quantum confinement 
effects, on surfaces. Of especial interest would be those systems in which these 
surface nanostructures form an ordered pattern. 

Among the technological issues behind surface nanostructuring, we can 
highlight: 

• New advances in magnetic storage technology [152]. 
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• The production of surfaces and thin films with well controlled shape at 
the nanoscalc for applications in catalysis [189]. 

• The design of nanoelectronic devices based on the use of quantum devices 
[3], for example for developing multispectral detector arrays by exploring 
novel detection techniques, nanopatterning and control of the production 
of self- assembled quantum structures. 

• The possibility of functionalization of these surface nanostructures that 
allows the selective attachment of specific molecules [92]. 

• The design of nanosensors based on surface-enhanced Raman scattering 
effects (SERS) or localized surface plasmon resonance [186]. 

• The use of these nanopattcrncd surfaces as templates to transfer these pat- 
terns to highly-functional surfaces [9] that can not be patterned directly. 
Also, the use of these patterns as templates can be applied to reduce fab- 
rication steps or to increase productivity. 

There are different approaches for surface nanostructuring. Among the 
so-called top-down methods we can mention lithographic techniques (nanoim- 
print lithography [76], nanosphere lithographic techniques [186], soft litho- 
graphic methods [184] and focused ion beam (FIB) techniques [117]). Other 
approaches are based on the use of SPM to induce nanostructures on a sur- 
face through different processes such as tip induced oxidation [18], tip induced 
e-beam lithography [155], dip-pen nanolithography [116], the application of 
strong electric fields between tip and sample [106], and the mere use of the 
probe as a stylus or pen to write at the nanoscale on the sample surface 
[156]. However, these methods present different limitations, such as proximity 
effects, low resolution or the need of parallel processing because of the (rel- 
atively) small processed area. In the case of SPM-based techniques, mainly 
Scanning Tunneling Microscopy (STM) and Atomic Force Microscopy (AFM), 
the use of an array of tips aims to compensate for this last limitation. 

Potential alternatives to overcome the limitations of the top-down ap- 
proach are provided by bottom-up methods, which are mainly based on self- 
organized processes. In this field, most efforts have been traditionally focused 
on the production of self-organized nanostructures occurring in semiconduc- 
tor heterostructure growth [159, 170] in which strain relief mechanisms take 
place. This simpler method seems to be a highly cost-efficient route towards 
large-scale arrays of nanostructures, although it presents some disadvantages 
such as enabling less control on structure size and shape than other methods, 
and the need to work under ultra-high vacuum conditions. 

In recent years, special interest is being paid to the study of self-organized 
nanopattern formation on surfaces by ion beam sputtering (IBS) techniques 
[55]. In general, two types of surface nanostructures can be induced by IBS: 
(a) nano-ripples and (b) nanodots. In both cases, the pattern formed by these 
nanostructures can have dimensions ranging from a few up to hundreds of 
nanometers. These patterns can be produced on different materials, amor- 
phous or crystalline, in just a few minutes and over areas of several square 
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millimeters. The diversity of materials processed and the similar morpholo- 
gies obtained indicate the universality of the process. In addition, it allows 
the eventual control of the induced nanostructures by changing the sputter- 
ing parameters, such as ion energy, dose, substrate temperature, ion incidence 
geometry, etc. Moreover, it can be used to produce functional surfaces and iso- 
lated structures [177]. Thus, IBS becomes a versatile, fast and cheap technique 
for surface nanopattering. 

The ability of IBS to induce submicro-structures on surfaces was reported 
more than forty years ago by Navez et al. [130]. In that work, they reported 
the production of periodic nanoripples on glass surfaces. In Fig. lb we show 
ripples obtained on a silicon surface immersed in an argon plasma. These rip- 
ples indeed remind us of those formed on sand dunes by the wind or on the 
sand bed close to the water edge by the water flow (Fig. la). As we will show 
below, this similarity is more than merely visual as the theories dealing with 
both types of (macro and nano) structures share many conceptual aspects. In 
the pioneering work of Navez and coworkers the authors did report important 
morphological observations, such as the change of the ripple orientation when 
the ion incidence polar angle is larger than a threshold value. Later, similar 
findings were reported for other materials, that have been reviewed elsewhere 
[21, 112, 127, 177]. Usually, the length scales involved in these patterns lie 
in the sub-micrometer scale, rather than in the nanometer scale. Moreover, 
the traditional field for ripple applications had been that of optical gratings 
[86]. From the theoretical point of view, the first advances for understand- 
ing IBS nanostructuring was made by Sigmund [154] as he showed that local 
surface minima should be eroded at a faster rate than local maxima (i.e. the 
sputtering rate depends on the local surface curvature), leading to a surface 
instability, which is the origin of the nanostructuring process. Based on this 
work, Bradley and Harper proposed later the first continuum model describ- 
ing ripple formation [13]. In recent years, other continuum models have been 
proposed [44, 111, 137]. These models account for different experimental be- 
haviors, such as the presence of absence of saturation for the ripple amplitude, 
ripple orientation, ripple dynamics as well as the existence or not of kinetic 
roughening. 

However, a clear turning-point occurred when in 1999 Fackso and cowork- 
ers [55] reported on the IBS production of GaSb nanodot patterns, which 
also display short-range order. This work shifted the focus of the research to 
the nanometer scale. The evident technological implications [34, 107] of the 
possibility to induce nanostructures on semiconductor surfaces on a relatively 
large area (several square millimeters) and in just a couple of minutes moti- 
vated further interest in these processes. Thus, nanodot production has been 
reported in different materials such as InP [64], Si [70], Ge [197], etc. With 
respect to the theoretical understanding of the nanodot IBS production, the 
theories must share most of the concepts with those developed for IBS induced 
ripples. An example is shown in Figures Ic and d, in which "dotted" struc- 
tures formed both on sand dunes and IBS nanodot structures are displayed. 
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Once more, this similarity appears that, rather than being just visual, is, in 
fact, deeper as will be showed in Sec. 4.3. 




Fig. 1. (a) Ripples on a sand dune in Morocco. Photograph courtesy of J. Rodriguez 
and E. Blesa. (b) 3.7 x 6.7 /im^ top view AFM image of a Si surface immersed in 
argon plasma, (c) "Dots" on a sand dune in New Mexico, USA. Copyright Bruce 
Molnia, Terra Photographies. Image Courtesy Earth Science World Image Bank 
http:/ /www. earthscienceworld.org/images. (d) 1x1 fim^ top view AFM image of a 
GaSb surface irradiated by 0.7 keV Ar^ ions under normal incidence. 

The production of dot patterns has changed the focus of the theoretical 
investigations. Now, for dot patterns the interest is more directed towards dot 
size and pattern order because of their technological implications. In the first 
case, it is interesting to know how the dot size changes with the different IBS 
parameters, specially ion energy, target temperature and sputtering time (i.e., 
fluence) . In the second case, it is necessary to know under which conditions the 
pattern order is enhanced. Thus, issues such as pattern wavelength coarsening 
and order enhancement have become relevant. Accordingly, new theoretical 
descriptions have been proposed [25, 35, 59, 87, 112, 123, 178]. 

The aim of this chapter is to give an overview of ripple and dot IBS 
nanopatterning that focus specially on those issues that remain open or, at 
least, ambiguous. This will be more evident for the case of nanodot patterns on 
amorphizable targets due, probably, to their relative novelty. In fact, we will 
treat specially the case of this type of patterns since previous reviews were 
devoted mainly either to ripple patterns [21, 112] or to the case of metals 
[177]. Stress will be paid to review the observed behaviors and the predictions 
of the main existing continuum models. When possible, we will contrast the 
experimental data with the theoretical predictions. 
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The chapter will be divided as follows. After this introduction and for the 
sake of completeness, we will present the basics of the Physics behind the ion 
sputtering process. Then, we will give an overview of experimental findings 
on surface nanopatterning by IBS on amorphous materials and single-crystal 
semiconductors. We will review briefly the results for the special case of single- 
crystal metal surfaces as they were already extensively reviewed by Valbusa 
and coworkers [177]. In the next section a review of the theoretical studies 
of these processes will be presented, with special emphasis on the various 
continuum models proposed so far. In the last section of the chapter we will 
propose possible applications of these nanopatterns and, finally, we will state 
the issues that still remain open, from our point of view, both theoretically 
and experimentally. 

2 Fundamentals of ion sputtering 

2.1 Introduction to ion sputtering 

Physical sputtering is defined as the removal of atoms from a solid surface 
due to energetic particle bombardment. This phenomenon was first reported 
by Grove [75] in 1852, although others had probably noted the effect while 
studying glow discharges. Grove studied the discharge generated by the tip 
of a wire positioned close to a polished silver surface and noted a ring-like 
deposit onto the silver surface when it was acting as an anode in the circuit. 
However, it was not until the early 1900s when the effect was attributed to the 
impingement of positive ions accelerated towards the cathode by the electrical 
field [160]. The sputtering process resembles the macroscopic sandblasting 
process but with ionized atoms instead of sand grains as bombarding species. 
In fact, the erosion efficiency with the incident angle of the abrasive particles 
[62] is analogous to that observed in ion sputtering experiments. 

The sputtering process is important both from a fundamental as well as a 
practical point of view. On the one hand, the imderstanding of the process is 
relevant to describe the basic interactions of ions with matter. On the other 
hand, as shown below, the process has found a broad range of applications. The 
improvement in experimental methods, as well as in the theoretical description 
of the process, has promoted the rapid maturation of the field during the last 
decades. 

2.2 Applications of ion sputtering 

The applications of ion sputtering to Surface Science and Technology are very 
diverse [127]. First, the removal of atoms from the surface can be used as 
an effective method for surface cleaning [167]. This may be aimed, for ex- 
ample, at removing the undesirable contamination layer present on a surface 
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(mostly containing oxygen and carbon impurities), which may affect the in- 
trinsic electrical or optical properties. This cleaning process may be crucial 
for the optimal performance of a semiconductor device [149]. 

The sputtered material from the target can also be transferred to a sub- 
strate for thin film, deposition (sputter deposition methods) [143]. This ap- 
plication of sputtering was first addressed by Pliicker in 1958 [138] after he 
observed the formation of a platinum film inside a discharge tube with a 
metallic mirror-like appearance. Sputter deposition is a non-equilibrium pro- 
cess where material ejection occurs at hyperthermal energies (in the few eV 
range) [165]. This energy range is much higher than in thermal evaporation 
methods, leading to coatings with better properties (higher density, etc.). The 
sputtering phenomenon is the base of many industrial coating activities, be- 
ing scalable from small pieces up to processing of large areas (it is commonly 
used to coat striictural window panels). It is relevant for the development of 
magnetron sputtering techniques [143], where high erosion/growth rates are 
obtained by confining the plasma discharge close to the cathode through a 
magnetic field. 

The ejected particles during sputtering can also be identified for compo- 
sitional analysis of the target material. This technique is known as secondary 
ion m,ass spectrometry (SIMS) [191] and is based on the principle that a small 
fraction of the sputtered particles are ionized. These particles are guided and 
detected by means of a mass spectrometer and provides very high depth res- 
olution (in the monolayer range) and elemental sensitivity down to parts per 
billion (ppb). 

The progressive etching of the target material through sputtering finds 
usage for sample thinning. Such a process is often incorporated into spec- 
imen preparation for transmission electron microscopy (TEM) where, after 
mechanical milling, grazing incidence ions impinge onto the remaining part of 
the sample until the desirable electron transparency is achieved (a few hun- 
dreds of nm for the energies commonly used for imaging). The progressive 
erosion of the surface is also an extended route for depth profiling with surface 
analytical techniques by a sequence of sputtering and analysis steps [133]. In 
this case, high depth-resolution can be achieved by erosion depths of a few 
monolayers but, as may occur in SIMS, by-products of ion bombardment such 
as atomic intermixing or preferential sputtering deteriorate the resolution of 
the analysis. 

The application of material etching through sputtering is more extended in 

the fabrication of integrated circuit devices. Here, pattern or feature formation 
taking part in implementing the design of the device is imprinted onto the 
substrate in several steps of ion beam etching or "drilling". This submicron 
architecture is defined by means of lithography, masks or FIB methods. 

Finally, as mentioned in the Introduction, sputtering can be used to alter 
the surface morphology due to its dependence on the local surface curvature 
[177]. In this case, depending on the conditions (angle of incidence, energy, 
etc.), it can induce either surface smoothening or roughening. These compet- 
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ing processes may be used to design the surface geometry required for certain 
applications. For example, smooth surfaces arc needed in optical components 
such as mirrors or lenses [67] or for tribological surfaces. For this purpose, 
sputtering should preferentially erode prominent topographies by, for exam- 
ple, grazing incidence bombardment, or promote ion-induced surface diffusion 
processes. On the contrary, roughening of a surface may enhance, for exam- 
ple, the surface activity (more effective surface is available), catalytic or other 
chemical processes, as well as control the biological response of the surface. 
Also, the competition between roughening and diffusion can lead to the evolu- 
tion of correlated features and induce the spontaneous formation of a regular 
pattern [177]. This self-organized pattern formation is the main subject of this 
chapter and will be described in detail in subsequent sections. 

2.3 Quantification of the sputtering process 

The sputtering process depends on several parameters such as the nature of 
the target, type of ion, ion energy, incidence angle, ion flux (rate of incoming 
ions towards the target per unit area) and total fluence (or dose) . For quantifi- 
cation of the ion-target interaction as a function of the different parameters, 
early investigations focused on the average number of target atoms ejected 
per incident ion. This magnitude is defined as the total sputtering yield (5) 
or, if measured in terms of a specific energy or angular interval, it is described 
as differential yield. However, there are other observablcs such as the effect 
on surface morphology or ion induced lattice damage that can be used for 
quantification of the sputtering effect. The understanding of the latter effects 
has only been possible with the recent development of proper characterization 
tools for surface imaging or structure determination. 

For multicomponent targets, the yield for each particular element should 
be taken into account. When the difference in 5" among species in the substrate 
is large, a preferential sputtering of the target elements with greater S occurs. 
This implies that the composition of the sputtered particles differs from that 
of the multi-component target and, as a consequence, induces an enrichment 
of the low- 5 component at the outermost layers of the target. 

2.4 Experimental measurements of the sputtering yield 

For a given ion and energy, there is no direct relation between the atomic num- 
ber of the target element and the resulting 5* [4]. However, there arc general 
experimental observations on the behavior of 5" as a function of other sput- 
tering parameters, mainly ion flux, energy and incidence angle. An important 
experimental observation is that ion erosion occurs even in the limit of low ion 
flux and that S is independent of the incoming flux. In addition, for a given 
projectile, S varies rather smoothly with incident ion energy until reaching a 
broad maximiun (in the 10-100 keV range) and, then, gradually dropping to 
zero for very high energies (in the MeV range) [190]. This behavior is found 
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for almost all projectile-target combinations. The variation of S with the ion 

incidence angle, 9, which is measured with respect to the surface normal, is 
considerable and increases from normal to oblique incidence until reaching a 
maximum around 60-70°, and then sharply dropping for glancing incidence 
angles [132]. The increase of S with 9 is explained by the higher energy de- 
posited in the surface region, a simple geometrical projection giving a cos~^ 9 
relation. The reduction of S at glancing incidence is due to the increase in 
the reflection of the incident ions [181]. The variation of the sputtering yield 
with the angle (relation between the maximum sputtering yield and the value 
at normal incidence) depends strongly on the ion/target combination, being 
larger for light ions. The general trends of as a function of ion energy and 
incidence angle are illustrated in Fig. 2. 




Ion energy (keV) Incidence angle, 6 (degrees) 



Fig. 2. (a) The observed energy dependence of the sputtering yield, S, of Si for 
normal incidence bombardment with different ions (adapted from Ref. [190]). (b) 
Normalized sputtering yield as a function of the polar angle of incidence for 1 keV 
Ar"*" on different materials (adapted from Ref [132]). 

It has been reported that S can considerably increase with ion fluence [4] . 
The steady state is reached only after the target has been eroded to a depth 
of the order of the projected ion range (maximum penetration depth of the 
ions in the direction normal to the surface). Obviously, this observation could 
be explained by the ion induced modification of the subsurface due to ion 
implantation. However, the analysis of the ejected particles also suggests that 
the yield increase may come from a mechanism of gas trap release [182]. 

Regarding other parameters, it is generally thought that the temperature, 
apart from small changes in the sublimation energy, does not influence sput- 
tering behavior significantly. Regarding crystallinity, the value of S may differ 
significantly as a function of energy among single crystal, polycrystalline or 
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amorphous targets of the same element. In single crystals, the yield may also 
depend on the surface orientation. A satisfactory explanation for these behav- 
iors relies on the concept of channeling (the collision probability is reduced 
if the ion is guided through a crystal channel), the contribution of this effect 
being less important in the low-energy range (below a few keV). 

In the previous description, energetic monoatomic particles impinging on 
the target surface have been considered. This is typically the case of sputter- 
ing experiments with noble gas ions. However, the sputtering process can be 
further complicated when the incoming ions react chemically with the target 
species or by bombardment with polyatomic or molecular ions. In the first 
case, chemical sputtering may also contribute to the increase of S by forma- 
tion of volatile compound products [145] . An example of this process has been 
observed by the increase of S during bombardment of Si with F or CI ions 
as compared with the values obtained for noble gas ions bombardment [166]. 
Conversely, when a non volatile compound forms by reaction of the incoming 
ion and the target element (for example, SiO formation by O bombardment of 
Si), the value of 5*18 generally reduced. In the case of polyatomic or molecular 
ions, it should be considered that the molecules fragment upon impact and S 
should account for the relative contribution of the individual constituents. In 
this case, a fair approximation considers that the individual ions have the same 
velocity as the incoming molecule and transfer their energy independently to 
the surface. 

Regarding the emitted particles, the spatial distribution of sputtered atoms 
is described generally as a cosine distribution. However, as a function of energy, 
this distribution can change into an under-cosine to over-cosine for low or high 
ion energy, respectively, or can be somehow more directional if the incident 
particle is inclined at a high angle with respect to the surface normal [134]. 
Usually, the direction of the sputtered cone is opposite to the incident direction 
and is thus called a forward-peaked distribution. 

Because of the energetic nature of the sputtering process, emitted atoms 
tend to have energies E well exceeding the thermal energy that would be 
consistent with the target temperature. The energy of the sputtered particles 
is typically peaked at a few eV and has a long tail that falls as E~'^ and goes 
up to the order of the incident ion energy [165]. The peak or average kinetic 
energy depend strongly on the mass of the target atom, as well as on the 
energy and mass of the incident particle. 

2.5 Theory of sputtering 

The many-body considerations of the sputtering process pose a formidable 
challenge to theoreticians. The first models to explain the sputtering process at 
atomic scale were based in terms of evaporation from a hot spot, i.e., "thermal 
spike" [160]. Although the angular distribution of sputtered atoms follows the 
distribution from a liquid phase and the presence of thermal spikes has been 
described by Molecular Dynamics (MD) simulations [52], this mechanism gives 
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very low S values due to the short-lived nature of the spikes [52] and it can 

not account for the flux independence in the sputtering yield. 

The theoretical understanding of ion sputtering processes was considerably 
improved by the introduction of collision cascades. This advance was achieved 
through a better description of the ion-solid interaction, which was triggered 
by the building of particle accelerators in the 1950-60s. When an energetic ion 
penetrates the surface of a solid target, it travels through the solid losing its 
energy through collisions with the nuclei and electrons of the target atoms. 
The range (maximum penetration depth) of the ions can be determined to a 
good approximation considering the nuclear {Sn = [dE/dx]n) and electronic 
{Se = [dE/dx]e) stopping independently [129]. Except for light ions (H, He), 
the nuclear stopping dominates at low energies (below a few keV) since the 
ions get closer to the atom. If the energy transfer to the target atoms is larger 
than the threshold given by the displacement energy, lattice displacements 
take place, which precisely occurs mostly when nuclear stopping dominates 
the energy losses. 

The concept of nuclear stopping was thus used by Sigmund [153] to de- 
velop a linear cascade model for sputtering. Here, the projectile shares its 
kinetic energy with target atoms initially at rest in a series of random binary 
collisions, a process in which fast recoils are created. The recoils knock on 
other target atoms creating a collision cascade in which more atoms are in 
motion, but progressively with slower speed. If the energy is deposited in the 
near surface region, atoms located close to the surface may gain sufficient 
energy and the appropriate momentum to leave the surface and be sputtered 
away. The collision cascade occurs in a time interval of picoseconds after the 
impact, when the recoil energies at the edges of the cascade have come below 
the threshold of displacement energy (of the order of a few tens of eV). The 
cascade is damped out by energy dissipation through several mechanisms, 
such as creation of phonons. 

By assuming a Gaussian energy distribution centered at the ion range, an 
analytical expression for the sputtering yield, S, can be deduced as follows 

S = aNSniE)/EB (1) 

Here, A'^ is the target atom density, Eb is the binding energy of a sur- 
face atom and a contains material and geometric parameters. The product 
NSn{E) represents the energy loss per unit length as the ion travels in the 
target. Sigmund's theory accounts satisfactorily for many of the observations 
in amorphous or polycrystalline elemental target materials, like the angular 
and energy dependences of S, as well as the angiilar and energy distributions 
of the sputtered particles. However, it is a continuum theory and, therefore, it 
does not address inherent inhomogeneity of materials, it can not be applied to 
crystalline structures since it neglects channeling effects and, finally, it does 
not consider other relevant issues such as the starting surface morphology. In 
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addition, purely repulsive forces between the ion and the target atoms are 
only valid for ion energies above a few keV. 

The correlation of sputtering with the surface morphology is a crucial issue 
for the content of this chapter. The effect of micro-rough surfaces was already 
posed by Sigmund in 1973 [154]. Differences in energy deposition, which are 
induced by topographical features, may cause a significant reduction of the 
local sputtering yield at local maxima and, inversely, an increase at local 
morphological minima. As a consequence, small irregularities on a relatively 
smooth surface may result enhanced by ion bombardment. This implies that 
microscopically flat surfaces are unstable under high-dose ion bombardment 
unless atom migration acts as a dominating smoothing effect. Conversely, 
sharp cones appear to be surprisingly stable under bombardment leading to 
surface nanostructuration by IBS. This prediction is the key-stone for the 
theoretical understanding of roughening or pattern formation induced by IBS, 
as will be shown in following sections. 

Regarding the fundamental description of the sputtering process, the 
shortcomings of the linear theory have been addressed by the development 
of computer simulations, providing a direct view of the collision events and 
incorporating the role of intrinsic material properties. In this case, discrete 
systems arc considered where the ions and target atoms interact in a deter- 
ministic way, like in MD simulations, or in a random nature, like in the case 
of binary collision approximation (BCA). In these models, the interaction is 
driven by repulsion forces that determine the ion trajectory in the solid. In 
the BCA method the interaction is treated sequentially until the energy of 
the projectile and recoils are thermalized. This method was implemented in 
the TRIM code [200, 201], where a screened Coulomb potential is considered 
for faster computation. The BCA is valid when the individual collisions are 
separated in space (fulfilled only for energy above a few keV) and fails when 
there is cascade overlap, as occurs for low ion energies. In this regime, attrac- 
tive forces among the particles become more relevant. Even at high energy, 
the BCA does not provide good description of the collision cascade as the 
particles slow down. 

In the low-energy range, a proper description should consider the simulta- 
neous movement of the projectile and recoils. Since typical processes during 
sputtering occur in the picoseconds regime, they can be properly described by 
MD simulations. In this way, experimental observations of crystalline targets 
have been successfully reproduced, such as the yield as a function of energy 
or the angular distribution of sputtered atoms [80]. However, MD calculations 
are time consuming in terms of computation and the interaction potential may 
not be accurate. In addition, the rather small length scales as well the short 
temporal window analyzed in MD simulations preclude a dynamical descrip- 
tion of the sputtering process by MD and make necessary the implementation 
of Monte Carlo simulations. 



12 Munoz-Garci'a, Vazquez, Cuerno, Sanchez-Garci'a, Castro, and Gago 



2.6 Experimental considerations for ion sputtering 

There are niany considerations that should be taken into account for precise 
and reproducible sputtering experiments. The large number of parameters 
involved and the variety of experimental configurations may imply a large 
dispersion of results that makes it difficult to present an unified picture. As 
an illustration of this complexity, most experiments for the determination of 
S before the 1970s are not reliable since they were not done under sufficient 
vacuum conditions [127]. Most experiments of surface nanostructuration by 
ion beam sputtering are done by means of low-energy broad ion beam guns 
(Kaufman, cold cathode, hollow cathode, etc.), medium-energy ions from ion 
implanters, focused ion beams (FIB) or directly by plasma immersion of the 
target. In most cases, ion beams are generated by a glow discharge and then 
extracted by means of proper beam optics. The nature of the discharge as 
well as the beam optics will define the range of ion flux and the ion energy 
that can be achieved. In all cases, several considerations should be taken into 
account such as energy spread, ion flux, vacuum environment or the eventual 
presence of multiple charge or mass ions. Moreover, a precise knowledge of the 
irradiation dose is required (direct ion current measurement with a Faraday 
Cup, collecting the sputtered material from the target or by measuring the 
etched volume by masking the target). Other relevant issues may be the unin- 
tentional temperature increase induced during the ion bombardment (target 
cooling is desirable) or the control over the starting material conditions (sur- 
face morphology, electrical properties, size, etc.) since small differences among 
experiments may lead to inconsistent results. 

Low-energy ion beams are difficult to handle in many cases due to space- 
charge blow-up. This is normally solved by neutralizing the beam with thermal 
electrons from a hot filament, by placing the ion source very close to the target 
or decelerating the energetic ions in front of the target. The presence of a 
ncutralizer may also be necessary for the processing of insulating materials. 
All these aspects may affect the beam divergence or beam focusing, which has 
been shown to play a relevant role in surface nanostructuring [43, 194]. 

Finally, broad ion beams are normally used to incorporate large-area pro- 
cessing. However, beam scanning in the case of small diameter beams, as in 
the case of FIB or ion implanter, can also rcsiflt in the formation of simi- 
lar patterns. In this case, the beam scanning may incorporate new effects in 
the nanostructuring process [38, 43], mainly driven by heating efi^ects on the 
target. 

3 Experimental observations of surface patterning by ion 
beam sputtering (IBS) 

As mentioned above, in the last years a relatively large amount of experimental 
works have appeared on the production of nanostructures and nanopatterns 
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by IBS. Roughly, the experimental reports can be divided into two groups: (a) 
those that merely report on the production of some IBS nanopattern, and (b) 
those works in which some systematic pattern analysis is done, related with 
the dependence on some of the experimental variables,. In this section we 
will try to summarize the most relevant results of the second type of works, 
with special emphasis on these open issues described in the Introduction. 
We will distinguish in this analysis two main categories depending on the 
nature of target material. Thus, we will describe separately the experimental 
findings for amorphous or amorphizable materials, such as single-crystalline 
semiconductors, and for crystalline materials that are not amorphized by the 
ion beam (i.e., metals). This division is grounded on the fact that continuum 
models are mainly based on the random distribution of the target atoms, 
which is the case for the first type of substrates. Although we will see that 
both types of materials share many characteristics when they arc subjected 
to IBS, the physical processes that actually take place can naturally lead to 
important differences in the pattern formation. An excellent recent review 
of IBS patterning on metal surfaces can be found in the previous work by 
Valbusa and coworkers [177], the experimental section will be mainly focused 
on IBS nanostructuring on amorphous and amorphizable materials. 

3.1 IBS patterning formation on amorphous or amorphizable 

surfaces 

Under this label, we include all the target materials that are amorphous (i.e., 
glass). In addition, we consider crystalline targets whose atom layers close to 
the surface become amorphous by the action of the impinging ions (i.e.. Si, 
GaSb, InP, etc). It is for these materials that Sigmund's theory [153] of ion 
beam sputtering is strictly applicable. 

In order to present the different results on the IBS nanostructuring of these 
materials, we will divide them into two main groups: (a) nanoripple patterns 
and (b) nanodot patterns. This is so because usually they are obtained under 
different experimental conditions and because the work on ripple patterning 
dates back to the 1960's whereas that performed on dot patterns is only 8 
years old. 

Ripple formation by off-normal ion incidence 

Ripple structures arc induced by IBS under off-normal conditions. However, 
this happens usually within a limited window of angles (for the Si case see 
[19, 183]). The ripple pattern has been reported in a huge amount of materials. 
To mention just a few examples, wc can quote observations on glass [130], 
Si02 [63, 114], Si [19, 49, 53, 71, 79, 176], Ge [28, 197], GaAs [48, 89, 110], 
InP [110, 51], highly oriented pyrolytic graphite (HOPG) [78], and diamond 
[47]. 
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Ripple morphology 

The ripple morphology is usually characterized by AFM, STM and TEM tech- 
niques because they provide enough resolution to reveal the ripple topogra- 
phy at the nanometer scale. When the ripple dimensions arc relatively large it 
can be assessed also by Scanning Electron Microscopy (SEM). Whereas from 
SPM techniques we can obtain morphological data such as ripple amplitude, 
surface roughness and even pattern ordering (see below), with TEM we can 
access structural properties such as the thickness of the ion-induced amor- 
phized layer as well any possible structure (for instance, bubbles) induced by 
the incorporation of the bombarding ions into the target [39]. One example 
of the information obtained by TEM is shown in Fig. 3. Here, we can observe 
that the amorphized layer thickness is quite different if the ripple slope faces 
the incoming beam or not, being larger for the face that is in front of the ion 
beam. Besides, some cavities corresponding to bubble-like features, due to the 
Ar+ incorporation to the target are clearly visible on the front slope. 



\ 




Fig. 3. (a) Cross-scctional TEM image of a Si surface sputtered witli 120 kcV Ar+ 
beam. The arrow indicates the ion beam direction. The horizontal bar corresponds 
to 500 nm. Micrograph courtesy of T K Chini. 

In contrast, the data typically obtained by SPM-tcchniqucs arc displayed 
in Fig. 4. In this figure we show a typical ripple pattern obtained on a silicon 
surface by Ar+ ion bombardment at ^ = 10°. In the cross-section we define 
two magnitudes that arc usually studied, namely, the typical wavelength. A, 
of the pattern and the amplitude of the ripple structures, A. In particular, 
the analysis of the former is very interesting as different models predict dif- 
ferent behaviors of A, for instance, with the ion flucncc. Thus, some models 
imply that A increases with sputtering time, that is, that the pattern coarsens 
with time. Therefore, the eventual coarsening of A is an object of study both 
experimentally and theoretically. 

Ripple orientation 

The dependence of the ripple orientation with 6 is one of the most relevant 
morphological properties of ripple formation by IBS. There is a critical angle, 
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Fig. 4. (a) 2 X 2 ^m^ AFM image of a Si surface bombarded by an Ar^ 1 keV 
ion beam at ^ = 10° for a fluence of 9 x 10^^ ions cm~^. (b) Cross section taken 
perpendicularly to the ripples in which A and A are indicated. 



6c, such that for 9 < 9c the ripples are perpendicular to the beam direction 
while for 9 > 9c the ripples run parallel to it. This fact, which was already 
observed on glass in the seminal work by Navez et al. [130] , is a classic behavior 
that was explained already by the theory proposed by Bradley and Harper 
[13]. This behavior has been observed for many target materials such as glass, 
Si02 [63] and HOPG [77]. One of these examples is shown in Fig. 5 for the 
case of fused silica. 




Fig. 5. (a) 1 X 1 fim^ AFM image of a fused silica surface bombarded by an Ar+ 
beam with 0.8 keV with $ = 400 pA cm~^: (a) for 20 minutes at 6* = 60° and (b) 
for 60 minutes at 6 — 80°. The dark arrow indicates the projection of the ion beam 
onto the target surface. Figure reprinted with permission from [63]. 



Ripple pattern: dependence on ion energy and type 

• Ion type — Ripple patterns have been produced by bombarding the target 
surface by different ions. The most frequently used species is Ar"*" due to its 
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low cost, inertness and relatively high mass. Among the noble gases, Kr+ 
and Xe"*" have also been employed [195]. Also, ripple patterns have been in- 
duced using beams of Cs+ [110], Ga+ [79], 0+ [108, 157] and N+ ions [157]. 
The two last cases imply, as explained before, that reactive sputtering ef- 
fects can take place. Therefore, ripple formation by IBS is quite general a 
process virtually independent of the target materials and bombarding ions. 

• Ion energy — The study of the dependence of the pattern wavelength with 
the ion energy, E, is quite interesting because it can be used to further 
check the consistency of the experimental erosion system with the assump- 
tions of Sigmund's theory [153] on which the different continuum models 
proposed so far are based. 

Regarding the energy of the ions, the following distinction is usually made: 
(a) low-energy range, which is normally applied to ion energies smaller 
than 2-3 keV, and (b) medium-energy range, which is applied to quite a 
wide range of energies ranging from 10 keV up to 100 keV or even higher 
values. Despite the wide range of energies, most of the works addressing 
the influence of the energy on the ripple pattern report, a qualitatively 
similar behavior is found, namely, that the typical wavelength increases 
with energy following a power dependence of the type A ~ E"^, where 
usually < m < 1. Figure 6 shows two examples of ripple production on 
silicon targets with low-energy (a) and medium-energy (b) Ar"*" ion beams. 
Clearly, the wavelength is quite different but the surface morphology is 
nonetheless similar. 




Fig. 6. (a) 3x3 /xm^ AFM image of a Si surface bombarded by a 1 keV Ar+ ion 
beam at 6* = 10° for ion dose of 2 x 10^^ ions cm~^. (b) 4 x 4 /xm^ AFM image of a 
Si surface bombarded by a 40 keV Ar^ ion beam at ^ = 50° for ion dose of 2 x 10^* 
ions cm~'^. The vertical bars indicate 150 nm and 900 nm, respectively. 



This behavior has been found for Si targets irradiated at low-energy for 
different ion species: [1, 176], Ar+ [195], and, Kr+ [195] and Xe+ [198], 
and also for sputtering experiments at medium-energy employing ions 
[91] and Ar+ ions under ion beam scanning conditions [38]. Also, this 
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dependence has been observed for Si02 targets bombarded by low-energy 

Ar"*" ions [175], for HOPG (bombarded by medium energy Ar"*" ions) [78] 
and diamond (bombarded by medium energy Ga+ ions) [47]. 
However, discrepancies arise when considering the values reported for the 
exponent m. It should be noted that this study becomes especially difficult 
for low-energy IBS experiments in which, usually, the sampled energy range 
is quite narrow as recently pointed out in [195]. Furthermore, sometimes 
it is difficult to assess whether the experimental data follow a linear or, 
rather, a different power law behavior. Thus, for different targets irradiated 
by low-energies ions, values of m in the 0.2-0.8 range have been reported 
[1, 175, 176, 195]. 

For the medium ion energy experiments, the energy range sampled is usu- 
ally wider due to the use of ion implanters. Therefore, the assessment 
of the quantitative dependence between A and E becomes more reliable. 
However, discrepancies still exist among the different values reported since 
m values in the 0.45-1 range have been obtained [38, 78, 91]. Particularly 
interesting is the case of HOPG targets: when they were bombarded by 
Ar+ ions a linear relationship was found, but when the ion was Xe"*" a 
power law behavior was observed with an exponent value rn « 0.7 [78]. 
Finally, quite a different behavior was found for Si targets bombarded by 
low energy Ar+ ions [16]. In this work A was found to decrease with en- 
ergy when the target surface was held at 717 °C, whereas no clear trend 
was observed for ripples produced at 657 °C. Also, an inverse relationship 
between A and ion energy was reported by Chini et al. for ion beam sputter- 
ing of Si surfaces without beam scanning [38] . Suppression or application 
of the ion beam scanning could influence the actual temperature at the 
target surface. This fact would be in agreement with a similar inverse be- 
havior observed for the experiments performed by Brown and Erlebacher 
at higher temperature [16]. 



Ripple pattern evolution with time or ion fluence 

The study of the ripple morphological evolution with irradiation time is usu- 
ally done through the dynamics of two magnitudes, the ripple wavelength and 
the siirface roughness. The former is related with the ripple lateral dimension 
and the latter with the vertical one (i.e., ripple amplitude). This study is quite 
important because eventual control of the pattern morphological properties 
would enable applications for technological purposes. 

• Ripple coarsening— The existence or not of ripple coarsening can be, as 
mentioned above, a touchstone for the continuum models. In fact, ripple 
coarsening is not predicted by the seminal BH theory [13] or some of its 
non linear extensions [44, 59, 111, 112, 178]. In contrast, it is predicted 
by more recent theories [25, 123]. Thus, the analysis of ripple coarsening 
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becomes a relevant issue. Moreover, since there are systems that display 
and others that do not show ripple coarsening, it is necessary to assess 
experimentally the differences between them in order to elucidate which 
physical phenomena are behind the coarsening process. This coarsening 
process typically reflects in a power law dependence such as A ~ where 
t is irradiation time (ffuence) and n is a coarsening exponent. 
An example of coarsening behavior (see Fig. 7) shows two AFM images of 
a Si(lOO) surface irradiated by 40 keV Ar+ ions at 50° for different times 
together with typical surface cross-sections of the surface morphology for 
both cases. In all images, the ripples run along the perpendicular direc- 
tion with respect to the projected ion beam direction (which runs along 
the horizontal axis of the top view AFM images) . From the AFM images 
coarsening of the typical ripple wavelength is already evident. This is better 
appreciated in panel (c), together with the clear surface roughening that 
also implies a clear increase in the amplitude of the ripple morphology. 




Fig. 7. Top view AFM images of a Si(lOO) surface irradiated by 40 IceV Ar"*" ions at 
50° for different times (ion fluences): (a) 20 x 20 /tim^, irradiation time 2 h (1.6 x 10^* 
ions cm~^) ; (b) 50 x 50 fira^, irradiation time 16 h (1.3 x 10^^ ions cm~'^). (c) Typical 
surface profiles taken along the projection of the ion beam of the samples irradiated 
for 2 h (bottom profile), 8 h (ion fluence of 6.5 x 10^*^ ions cm~^; middle profile) and 
16 h (top profile). 



For the case of Si targets, we find ripple coarsening for low-energy ion irra- 
diation experiments with 0.5 keV Ar+ ions impinging at 60° [17] . Moreover, 
this coarsening process has been observed for a relatively wide substrate 
temperature range, from 600 °C up to 748 ° C [16]. Also, a coarsening 
process was observed for 1.5 keV Ar+ ions impinging at 45° and room 
temperature [103] and 1 keV ions impinging at 52° [108] for both high 
and low ion fluxes. 

In contrast, coarsening was not observed for irradiation experiments per- 
formed at 582°C and at 67.5° with 0.75 keV Ar+ ions [53]. This was also 
the case for experiments done by bombarding the Si surface at ^ = 15° 
with a 1.2 keV Ar+ ion beam. The same behavior was found when the ion 
species were Xe+ or Kr+ [195]. 
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For the medium-energy ion range, ripple coarsening has been reported for 

Si surfaces bombarded at 6* = 45° by 20 kcV and 40 keV Xe"*" ions where 
target temperatures were maintained between 100-300 K [19]. Similarly, it 
has been observed for Si targets irradiated either by 60 keV Ar+ ions at 

e = 60° [49] or by 30 kcV Ga+ ions at 61 = 30° [79]. 

In other materials, ripple coarsening has been observed for HOPG and 
diamond surfaces. In the first case, the irradiation process was performed 
by 5 kcV Xc+ ions at = 60° and 9 = 70° [77] whereas in the last case 
Ga+ ions impinged at 9 = 60° with energies of 50 keV and 10 keV [47]. 
Also, ripple coarsening has been reported for fused silica bombarded by an 
0.8 kcV Ar+ ion beam at 9 = 60° [63] and for glass targets irradiated by 
0.8 keV Ar+ ions, which were generated in a defocused electron cyclotron 
resonance plasma, with an angle of incidence of 35° [174]. Finally, ripple 
coarsening was also observed for InP targets irradiated at 9 = 41° by 0.5 
keV Ar+ ions [51] but not for GaAs surfaces bombarded at 9 = 41° by 
10.5 keV Ot ions [89]. 

With regard to the value of the exponent n different values have been 
reported. Thus, on silicon targets a value of 0.5 for irradiation at 30 keV 
has been reported [79]. Also for high energies, 60 keV, two regimes with 
values ni = 0.64 and ni = 0.22 were observed. In contrast, for low-energy 
irradiation experiments an exponential dependence, rather than a power- 
law, was found for a relatively wide range of substrate temperatures [16]. 
When 0.8 keV Ar+ ions were employed to irradiate fused silica [63] and 
glass surfaces [174] coarsening exponents of 0.15 and 0.95 have been found, 
respectively. 

Finally, only for 60 kcV Ar+ ion irradiation at 9 ~ 60° of GaAs surfaces, 
a disordered dot morphology has been reported previous to, and later co- 
existing with, the ripple morphology [48]. For long sputtering times (i.e., 
ion fluences of 3 x 10^* ions cm^^) only the ripple morphology remained, 
without any nanodot structure superimposed. 

• Surface roughening — The study of the surface roughness, W defined 
as the mean square deviation of the local height with respect to its mean 
value — of the patterns can be very useful for both technological (e.g. de- 
veloping metal surfaces for SERS applications) and fundamental purposes. 
In principle, W should be proportional to the ripple amplitude (i.e., the 
peak to valley height difference). A, in case the patterns were perfectly 
periodic. However, in real patterns there arc height fluctuations among 
ripples that imply that A and W are not completely equivalent. Although 
most part of the studies deal with W, some of them analyze A rather than 
W. 

The most frequently observed behavior is that W initially increases steeply 
(usually increasing exponentially with ion fluence or sputtering time) to 

either saturate or grow at a slower pace, usually following a power law 
dependence such asW^t'^.ln. the latter case, it is interesting to measure 
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the value of f3 since it can be contrasted with predictions from theoretical 
models. 

For the first case, i.e., exponential increase followed by saturation, we can 
mention experiments on Si surfaces irradiated by low-energy Ar+ ions [53] , 



The second case, i.e., sharp increase followed by power law behavior, has 
also been observed in many systems: Si targets bombarded by medium- 
energy Ar+ ions [49, 91] and also for HOPG surfaces bombarded under 
similar conditions [78]. For Si targets irradiated by 60 keV Ar+ ions an 
initial value of /3i = 0.76 was reported, although this relatively high value 
could be compatible with an exponential increase, whereas for longer times 
P2 = 0.27 was found [49]. In addition, for 16.7 keV ions /3 = 0.38 was 
observed after the initial sharp increase ofW [91]. Also, for HOPG surfaces 
irradiated by 5 keV Xe+ ions a (3 value compatible with the Kardar-Parisi- 
Zhang [88] universality class was reported [77, 78]. 

Finally, there are works where only power-law behaviors were reported. 

Most of these systems studied present a /3 value such that 0.45 < /? < 1 
[19, 49, 63, 174] (in the second case, for low ion flux conditions). These 
behaviors could be due to partial analysis of the initial exponential in- 
crease that, when analyzed in a limited temporal range, can be analyzed 
in terms of a power law dependence with an exponent /? > 0.5 as remarked. 

• Shadowing effects An important issue that should be taken into ac- 
count for studying ripple evolution is geometrical shadowing. These ef- 
fects appear for long sputtering times and when relatively large 6 values 
are employed. The important role of shadowing effects was highlighted in 
[19]; subsequently. Carter gave a simple estimation of the conditions under 
which these effects begin to operate [20]. In particular, he proposed that 
shadowing operate for ripples with an amplitude (here taken as propor- 
tional to W) to wavelength ratio such that 



In equation (2) as 9 increases (i.e., approaching grazing incidence) the right 
hand side making shadowing effects more likely to appear. These effects 
have been considered recently for Si irradiation by 60 keV Ar+ ions at 
6 — 60° [49, 50]. According to (2), shadowing effects should appear when 
W/X > 0.09, which occurs after 800 s of irradiation under the sputtering 
conditions described in [49]. Interestingly, this value is very close to the 
threshold time for which ripple coarsening begins to be observed. Thus, 
shadowing effects can influence largely the ripple dynamics. This fact can 
be important because usually shadowing is not incorporated into contin- 
uum models. Moreover, the latter are usually derived under a small slope 
approximation and the abrupt morphologies generated when shadowing 
processes appear cannot be described under such approximation. A re- 
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W^/A>27rtan(7r-6l). 



(2) 
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markable exception is the continuum equation proposed in [35] that seems 
to correctly describe steep surface features. 



Ripple pattern dependence on target temperature 

The study of the pattern evolution with target temperature can contribute 
to further increase our knowledge of the main physical mechanisms governing 
the IBS pattern formation. In particular, temperature can affect the surface 
diffusivity, which can lead to changes in the morphology of the pattern. How- 
ever, studies are scarce due, probably, to their experimental complexity. For 
low-energy ions there arc two reports on silicon surfaces irradiated by Ar^". 
In the first one, in which 0.75 keV Ar"*" ions impinged aX 9 = 67.5° onto the 
surface [30, 53], A coarsened with the target temperature in the 460-600° C 
range following an Arrhenius law [112] 

^~5n7^exp(^^j. (3) 

These data are shown in Fig. 8. The value of A obtained at the lowest 
temperature was measured by AFM because it was out of range of the light 
spectroscopy measurements. This analysis led to a value for the activation 
energy, {AE = 1.2 ± 0.1 eV) for surface mass transport on ion-bombarded 
Si(OOl). 



m 




Tempemtine ( °C ) 



Fig. 8. Ripple wavelength A versus 

target temperature for a Si(OOl) sub- 
strate irradiated at 6* = 67.5° by 
0.75 keV Ar+ ions. Open squares cor- 
respond to light scattering measure- 
ments, whereas the filled circle was ob- 
tained by AFM. The solid line corre- 
sponds to a fit following the Arrhe- 
nius law (3). Figure reprinted from [53] 
(http: / /link. aps.org/abstract /PRL /v82 
/p2330) with permission 



The same behavior was found for Si(lll) irradiated at 60° by 0.5 keV 
Ar+ ions in the 500-750° C range [16]. This study led to AE = 1.7 ± 0.1 
eV. Similarly, for Si02 targets irradiated by 0.5-2 keV Ar+ ions an Arrhenius 
behavior was observed for T > 200° C [175]. 



22 Munoz-Garci'a, Vazquez, Cuerno, Sanchez-Garci'a, Castro, and Gago 



Finally, there are two reports on the variation of the pattern with target 
temperature for medium-energy IBS experiments. In the first one, HOPG 
surfaces were irradiated by 5 keV Xe"*" ions in the 573-773° C range [78]. 
For this system the ripple wavelength followed also the Arrhenius law giving 
a value of AE = 0.14 cV. The second experiment was on bombardment of 
GaAs targets by 17.5 keV Cs"*" ions at = 25° [110], in which the ripple 
wavelength was analyzed in the — 100° C range. The data obtained have been 
later analyzed in [113]. For T > 60°, the behavior was well described by Eq. 
(3) with an activation energy for surface self-diffusion of AE = 0.26 eV. For 
lower temperatures, a slight decrease of A was observed following a T~^l'^ law. 
However, within error bars, these experimental data are also compatible with 
a temperature independent behavior, in agreement with models for effective 
surface diffusion effects of erosive origin [111, 112]. In any case, the present 
experiment [110, 113] provides a clear example of the existence of different 
temperature regimes for ripple formation under IBS. 

In contrast, Carter and Vishnyakov did not find any change of A with 
temperature in the sampled 100-300 K range when Si targets were irradiated 
by 10-40 keV [19]. These findings are consistent with the existing theories since 
for relatively low temperatures ion-induced surface diffusion processes, which 
not depend on the target temperature, dominate over thermally activated ones 
[111]. 

Ripple pattern dependence on ion flux 

The ripple pattern morphology, in particular its wavelength, can depend also 
on the ion flux, ^, i.e., the number of incoming ions per area and time units. 

Among the scarce studies of this behavior, most of them did not find any 
change of the ripple wavelength with ion flux. This was the case for Si surfaces 
bombarded by either 1.5 keV ions at 6* = 40° [176] or by low-energy Ar+ 
ions at different ion fluences [16] or by 2 keV Xe+ ions at ^ = 20° [198]. In 
addition, the same behavior was found for fused silica targets irradiated at 
9 = 60° at different angles by a low-energy Ar+ ion beam [63]. Also, no change 
of A with the ion flux was obtained for diamond surfaces irradiated by 50 keV 
Ga+ FIB a,t 6 = 57° [47]. All these studies were done at room temperature, 
except for the one on fused silica that was performed at 12° C. 

In contrast, for Si surfaces bombarded by 0.75 keV Ar+ ions at 67.5 degrees 
to normal in the 500-600° C range [53], A was found to decrease with as 
\ Also, a decrease of the ripple wavelength with ion flux was reported 

for Si surfaces bombarded at room temperature by 1 keV ions at ^ = 52°, 
although in this case the quantitative dependence was not addressed [108]. 

Ripple pattern order 

The size of ordered ripple domains is another essential property of these pat- 
terns for potential technological applications. However, it is somehow difficult 
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to assess. In principle, there may be several methods to evaluate it. A first 

one is based in the data obtained by AFM. From these data, it is straightfor- 
ward to obtain the Power Spectral Density (PSD) of the surface morphology, 
h{-r, t). The PSD is defined as P_SD(fc, t) =_ (/i(k, t)/i(-k, t)), where /?,(k, t) is 
the Fourier transform of /i(r, t) — h{t) with h{t) being the space average of the 
height and k = |k|. This PSD curve usually presents a peak denoting the exis- 
tence of a characteristic mode whose associated length scale is identified with 
the ripple wavelength A. It has been proposed that the pattern lateral corre- 
lation length, Q, which gives an estimation of the average size of the ordered 
domains, can be obtained from the full width at half maximum of the PSD 
peak [193]. This method was employed by Zibcri et al. [195] for estimating 
the range of order of ripple patterns produced by 1.2 keV Ar+, Kr+ and Xe"*" 
ions at ^ = 15° on silicon surfaces. Although in these experiments coarsening 
was not observed, (thus A was constant for all the ion fluences), in all cases 
C was observed to increase with the ion fluence (especially, for the largest ion 
fluences C > llA was found). In this work, the authors also studied the change 
of C with the ion energy for the three different ions employed. Whereas for 
Xe"*" and Kr"*" ions C, and A were found to increase in the same way with ion 
energy, for the case of Ar+ ions a maximum for the C/A ratio was observed 
for an ion energy of 1.2 keV. It should be noted, that for the largest energy 
employed in all cases, 2 keV, C ~ lOA, irrespectively of the ion species. 

The previous method for assessing the pattern order degree has the dis- 
advantage of the local character of SPM techniques. However, there is an- 
other method, based on grazing incidence diffraction (GID) or small angle 
scattering (GISAXS) synchrotron techniques, that provides better sampling 
statistics. However, to our knowledge, there is no published work using this 
technique for such purposes on amorphous materials where it has been used 
ripple crystallinity assessment [82]. 

Ripple propagation 

The transversal collective motion of ripples has only been studied in two works 
by simultaneous real time monitoring of the ion-induced ripple morphology 
by SEM. In the first one [79] a Si surface was irradiated by 30 keV Ga+ ions 
and 6 = 30°, it was found that ripples initially propagated along the ion beam 
projection direction with a velocity of v = 0.33 nm s^^ to slow down later on 
as they coarsened. In the second case, in which glass surfaces were irradiated 
by 30 keV Ga"*" ions [2], ripples also propagated along the projection of the 
ion beam direction. In this case, the ripple evolution was shown in real time. 
Moreover, ripples did not initiate the propagation until most of them were 
completely formed, to finally reach an uniform propagation velocity. 
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Nanodot patterning in amorphous/amorphizable materials 

Dot IBS nanopatterns are produced when the anisotropy caused by the obhque 
incidence of the ion beam is suppressed. Basically, there are two ways to elim- 
inate this anisotropy: (a) by IBS under normal incidence, which is the most 
frequently employed technique [55, 70]; (b) by IBS under oblique incidence but 
with simultaneous rotation of the target [64]. More recently, nanodot patterns 
have been produced on Ge surfaces by a 2 keV Xe"*" ion beam impinging on 
the Ge surface at 6* = 20° [197], possibly related with the role of the critical 
angle 6 = 6c mentioned at the beginning of Sec. 3.1. 

These patterns arc usiially characterized by a highly uniform dot size dis- 
tribution and short-range in-plane ordering. These two properties make them 
very interesting for potential technological applications. Although the first 
report on the production of such IBS nanopatterns [55] is relatively recent, 
dating from 1999, many groups are investigating the mechanisms leading to 
their formation. Thus, up to now, these patterns have been produced in dif- 
ferent materials: GaSb [12, 55, 66, 185], InP [64], InAs [68], InSb [57], Si 
[70, 71, 196], and Ge [197]. 

As occurred for the ripple patterns, the formation of the dot patterns 
under different experimental conditions and for a relatively wide range of 
materials suggests that this process does not depend on the specific ion-target 
interactions. 

Another important issue regarding the target material is its crystallinity. 
Thus, it has been proved that, for GaSb targets, IBS nanodot patterns can be 
produced both on crystalline [55] and amorphous surfaces [58]. In addition, 
these patterns have been produced by IBS on both Si(lOO) and Si(lll) surfaces 
[72]. In this work, it was found that, although the surface crystallinity does 
not affect the pattern formation, it can affect to some extent the pattern 
dynamics. In particular, it was observed that Si(lll) surfaces have a faster 
dynamics in terms of pattern coarsening (see Fig. 9) and ordering than Si(lOO) 
surfaces. Again, the fact that the surface crystallinity does not determine the 
formation of the pattern is compatible with Sigmimd's theory [153] for which 
the surface is considered as amorphous, cither originally or as induced by the 
ion beam action. 

In the following, we will rcvi(nv the main experimental findings on the 
properties of these patterns depending on the different experimental parame- 
ters. One of the most studied issues is the variation with physical parameters 
of the basic pattern length scale, A, which now corresponds to the average 
dot-to-dot distance, usually proportional to the dot size. It is worth noting 
that the interest in these studies on the dot shape and size lies in the re- 
quired control of these parameters, particularly the dot size, for developing 
technological applications. 

It should be noted that, whereas the dot size is usually affected by the 
AFM tip (AFM being the routine technique for characterizing the surface 
morphology), the dot-to-dot distance is not usually affected by tip convolution 
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Fig. 9. Dependence of A with sputter- 
ing time obtained by GID for Si(OOl) 
(•) and Si(lll) (o) surfaces irradiated 
under normal incidence by 1.2 keV Ar"*" 
ions. Note the faster coarsening dy- 
namics for the IBS pattern induced on 
Si(lll) surface. 



effects [65]. The pattern wavelength is usually determined from the radially- 
averaged PSD of AFM images as is shown in Fig. 10. In panel (a) of this figure 
we display a typical AFM image of a nanodot pattern induced onto a Si(lOO) 
surface. Panel (c) shows the radially-averaged PSD function corresponding 
to this image, in which we can observe a dominant peak corresponding to 
the basic wavelength A. At higher k values we can detect other minor peaks 
or shoulders indicating the high lateral ordering and size homogeneity of dots 
[199]. Conversely, panel (c) also illustrates the power-law behavior of the PSD 
at small fc, which signals height disorder between dots at long distances. This 
behavior is very frequently found for this type of experiments. 




Fig. 10. (a) 1 X 1 /xm^ AFM image of a Si(OOl) substrate irradiated under normal 
incidence by 1.2 keV Ar"*" ions, (b) Two-dimensional auto-correlation function in 

which the short-range hexagonal order is observed. The horizontal bar corresponds 
to 68 nm. (c) Radially-averaged PSD for image of (a). The main peak corresponds 
to the k = 2'n/X value corresponding to the dominant pattern wavelength. 



Similarly to the case of ripples, TEM analysis can provide us useful in- 
formation regarding the morphology and structure of the dot patterns. Thus, 
different dot morphologies obtained by TEM are presented in Fig. 11. In (a) 
we observe the conical morphology of crystalline GaSb dots produced under 
normal irradiation and target rotation [69]. In contrast, in (b) we observe 
GaSb dots with a sinusoidal shape obtained under irradiation at ^ = 75° and 
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target rotation. In these cases, the amorphous layer thickness was ~ 4 nm 
[69]. Finally, in (c) we observe Si nanodots produced under normal irradi- 
ation without target rotation displaying, rather, a lenticular shape with an 
amorphous layer thickness of cr; 2 nm [70] . 





Fig. 11. Cross-sectional TEM im- 
ages of (a) GaSb dots nanopatterns 
obtained under 0.5 keV Ar"*" irradia- 
tion with a fluence of 10^^ ions cm~^ 
at normal incidence with target rota- 
tion; (b) GaSb dots nanopatterns ob- 
tained under 1.2 keV Ar"*" irradiation 
with a fluence of 10^^ ions cm~^ at 
6 = 75° with target rotation; (c) Si 
dot nanopatterns obtained under 1.2 
keV Ar"*" irradiation with a fluence of 
9 X 10^^ ions cm~^ at normal inci- 
dence. The horizontal bars correspond 
to 50 (a), 25 (b) and 30 nm (c), respec- 
tively. Figures (a) y (b) taken from 
[68] with permission. 



Nanodot pattern dependence on ion energy and type 

• Ion type — To date, all the experiments in which nanodot patterns have 
been reported are done with low-energy Ar+ ion beams, except for one 
work in which Ne"*", Kr"*" and Xe"'" ion beams were also employed [199]. 
In this work oblique ion beam incidence and simultaneous target rotation 
were employed. Although some characteristics of the pattern differed, the 
authors concluded that there was almost no difference in the morphologi- 
cal evolution of the mean size of dots when using different ions. In contrast 
to the experimental findings of the same group on ripple formation, the 
use of Ne"*" ions did lead to the production of nanodot patterns with an 
experimental behavior similar to that observed when Ar"*" ions were em- 
ployed. Recently, experiments have been performed using 1 keV ions 
[168]. However, in this work a sort of nanodot chains was produced that 
seem more similar to nanoripples with a superimposed nanodot morphol- 
ogy- 

• Ion energy — Regarding the energy range used in the nanodot experi- 
ments, there is a marked difference with respect to the studies realized on 
nanoripple IBS production (see above). Namely, for nanodot experiments 
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only low-energy ions have been used, with E <2 keV. As remarked earlier 
in Sec. 3.1, working in this energy range can lead to relatively large errors 
in the determination of the power-law dependence A ~ , especially in 
the value of m [66] . However, this was not the case for the study by Facsko 
ct al. for IBS nanodot patterns on GaSb since they sampled 15 energy 
values within this range [57]. They obtained a value of m = 0.5 ± 0.02. In 
the same study the authors estimated, a similar m value for InSb surfaces 
from three experimental data points. 

In Fig. 12 we present our results on GaSb surfaces irradiated at different 
ion energies. Clearly, the dot size coarsens with ion energy. In fact, our 
experimental data arc consistent with an exponent of m = 0.5 as in the 
case of Facsko and coworkers [57]. Also we must note the evident surface 
roughening with ion energy given that the vertical scale is the same in all 
three AFM images. 

Another experimental study was done by the group of Frost and coworkers 
for Si surfaces with Ne+, Ar+, Kr+ and Xe"*" oblique ion beams onto rotat- 
ing Si targets [199]. In all cases, they found an increase of A with ion energy. 
Unfortunately, they did not estimate the value of m. However, from a vi- 
sual inspection of Fig. 6 of [199] it is clear than rriNe ~ "^Ar > ^riKr > ^^^Xe• 




Fig. 12. (a) 1 X 1 fiui^ AFM images of a GaSb substrate irradiated with Ar ions 
under normal incidence with an approximate fluence of 10^*^ ions cm~'^ at: (a) 0.3 
keV; (b) 0.7 keV; (c) 1.2 keV. The vertical bars correspond to 200 nm. Note the 
surface roughening and dot-size coarsening with increasing ion energy. 
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Nanodot pattern evolution with sputtering time or ion fluence 

The study of the influence of sputtering time (i.e., ion fluence) on the pattern 
morphology has attracted the interest of the researchers already since the 
seminal work by the group of Facsko and coworkers [55]. As occurred for 
the ripple morphologies, two main dynamics are studied: that of the pattern 
wavelength and the evolution of surface roughness. 

• Dot nanopattern coarsening — Basically, the same behavior is observed for 
many of the various experimental systems: initially, A increases to satu- 
rate afterwards. However, the dynamics of this process is quite different 
depending on the target material and ion current density. Thus, for GaSb 
surfaces [12, 185] and InP [64, 66] the saturation regime is attained for 
ion dose close to 10^* ions cm~^. It should be noted that this saturation 
was attained for an ion dose of 1.7 x 10^^ ions cm~^ when InP targets 
were irradiated under normal incidence by Ar"*" ions without rotation, but 
with a flux six times smaller [169] than that employed in [64]. For Si sur- 
faces saturation takes place for a considerably larger ion dose, at 4 x 10^^ 
ions cm^^ [72]. In another study where the Si surface was intentionally 
seeded with molybdenum, A can be estimated to saturate at an ion dose 
of 7 X lO^'' [135]. The different dynamics of GaSb and Si surfaces under 
normal incidence Ar"*" irradiation and rotating InP targets under oblique 
Ar"*" bombardment is shown in Fig. 13, where results obtained on these 
systems are displayed. 
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Fig. 13. Dependence of the charac- 
teristic pattern wavelength with ion 
fluence for GaSb, $ — 0.8 rriA crn~^ 
(□), InP at ^ = 0.15 mA crn"^ (■) 
and Si at * = 0.24 mA cm"^ (•) 
surfaces. The data for GaSb and InP 
have been adapted from [12] and [64], 
respectively. 



Mcasiirement of the exponent n in the power law dependence A ^ before 
saturation has been done for four systems: (a) for Si surfaces irradiated 
under normal incidence by 1.2 keV Ar+ ions with = 0.24 mA cm~^ 
Gago et al. foimd n ^ 0.2 [70]. (b) For GaSb targets sputtered by 0.5 
keV Ar"'" ions under normal incidence with = 0.8 mA cm~^ Xu et al. 
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reported n = 0.14 ± 0.03 [185]. (c) For rotating InP targets irradiated by 
oblique 0.5 kcV Ar+ ions with = 0.15 niA cm^^, Frost ct al. reported 
a value of n = 0.26 ± 0.04 [64, 66]. (d) For InP targets irradiated under 
normal incidence by 1 keV Ar+ ions with = 0.0233 mA cm~^ , Tan et 
al. reported a value of n = 0.23 ± 0.01 [169]. 

In addition, for Si surfaces irradiated by 0.5 keV Ar"*" ions Ludwig et al. 
did observe a coarsening process with sputtering time but they did not es- 
timate the value of the coarsening exponent [109]. However, they did not 
observe the saturation regime up to ion doses of 4.8 x 10^^, which agrees 
with other experimental reports on Si surfaces. Besides, the authors did 
find that the coarsening exponent value increased with the ion energy in 
the 100-200 eV range. However, the opposite behavior, i.e., absence of 
coarsening, was observed by the group of Frost for rotating Si targets ir- 
radiated by oblique beams of Ne+, Ar+, Kr+ and Xe+ ions [196, 199]. 

• Surface roughening — The surface roughness gives us a measure of the 

nanodot height as well as of the height fluctuations among dots. For GaSb 
surfaces [12], W initially increases rapidly with sputtering time to reach 
later a maximum value. In this first region Xu and Teichert obtained /? = 
0.87 ± 0.12 [185]. Again, this high value could be an indication that the 
surface roughness increases exponentially rather than follow a power law. 
For longer times the surface roughness decreases to attain a saturation 
value. 

For the case of InP targets irradiated under normal incidence conditions a 
first regime for which /3i = 0.74 ± 0.03 was reported [169]. This regime led 
to another one in which /32 = 0.09 ± 0.03, also roughly compatible with 
roughness saturation. In contrast, for rotating InP targets [64], the surface 
roughness increases for the whole temporal range that was sampled. A first 
value /?! = 0.8 ± 0.1 was measured, once more possibly compatible with 
an exponential time dependence, while for longer times P2 = 0.27 ± 0.06 
was obtained. 

For Si targets, both fixed and rotating, a similar behavior of the roughness 
has been reported, namely a sharp initial increase followed by saturation 
[70, 135, 199]. In the latter study, the behavior observed for Ar+, Kr+ and 
Xe"*" ions was found to be consistent with an initial regime during which 
the roughness increased exponentially. 

The main findings of the above works are displayed in Fig. 14 where we 

plot the evolution of W with sputtering time for the GaSb, InP and Si 
systems. In this plot we have represented the x-axis in logarithmic scale in 
order to display the three systems in a single graph, which have different 
dynamics. It can be appreciated that, although the temporal evolution is 
different for each system, there is always a sharp initial increase of the 
roughness before reaching either a stationary value or a regime with a 
slower growth. In the inset we display the same plot only for the initial 
stages of the sputtering process (i.e., t < 170 s); now the y-axis is the 
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one in logarithmic scale so that for all three systems the initial roughness 
seems to increase exponentially with time as the straight lines suggest in 
the plot . 
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Fig. 14. Surface roughness vs sput- 
tering time for IBS dot patterning 

of GaSb irradiated under normal in- 
cidence (■), InP bombarded under 
oblique incidence and simultaneous 
target rotation (o) and Si irradiated 
under normal incidence (•). Note that 
the temporal axis is in logarithmic 
scale. Inset: same graph but restricted 
to short sputtering times, i.e., t < 170 
s. Note that now the y-axis is in log- 
arithmic scale. The solid straight lines 
are guides to eye to indicate the expo- 
nential dependence of the initial rough- 
ness increase for the GaSb, InP and Si 
systems, respectively. 



Finally, we have also studied [70] , how the surface roughness changes with 
the substrate temperature for Si targets irradiated under normal incidence. 

We observed that W, which was constant for temperatures up to 400 K, 
decreased to reach a saturation value at 550 K where the pattern vanished 
[73]. 



Nanodot pattern dependence on target tem.perature 

To the best of our knowledge, there are four studies on the dependence of 
the pattern morphology on the target temperature, each one on a different 
material. Thus, we will describe separately the main findings of these studies. 

1. GaSb surfaces irradiated under normal incidence by 0.5 keV Ar+ ions [57]: 

for this system the pattern wavelength did not change with target tem- 
perature between —60° C and 60° C. This behavior was interpreted as a 
confirmation that the main relevant smoothing process is nonthermal un- 
der these experimental conditions; the main relax;ation process is, rather, 
due to ion- induced effects [111]. 

2. Rotating InP surfaces irradiated by 0.5 keV Ar+ ions impinging at 30° 
[66]: in this case, quite a striking complex behavior was found since pattern 
symmetry changed in the temperature range between 268 K and 335 K 
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from short-range hexagonal to square patterns symmetry. In addition, the 
characteristic wavelength increased with temperature [68] . This behavior 
is shown in Fig. 15. 




Fig. 15. AFM images of a rotating InP target irradiated at ^ = 30° by 0.5 keV 
Ar+ at ^ = 300 /xA cm~^ at different target temperatures: (a) 268 K, (b) 335 K. 
The corresponding two-dimensional autocorrelation functions are displayed in (c) 
and (d), respectively. Figure reprinted from [68] with permission. 

3. For Si targets irradiated by 1.2 keV Ar+ ions [73] was observed, by both 
AFM and GISAXS, the nanopattern wavelength to be a constant up to 425 
K, and then to decrease in the 425-525 K range. For higher temperatures 
the pattern vanished and the surface became featureless. This behavior is 
not explained by any of the existing theories on IBS nanostructuring. In 
Fig. 16 we display AFM images of IBS induced dot patterns at 300 K (a) 
and 425 K (b) showing qualitatively how both A and the dot size become 
smaller with increasing substrate temperature. This shrinking process be- 
comes evident in panels (c) and (d) in which we show the PSD and GID 
curves measured on both patterns, respectively. In both cases, the main 
peak shifts to higher k and q values, i.e., the characteristic length scale 
diminishes as substrate temperature increases [73]. 

4. Fixed InP surfaces irradiated under normal incidence [169] by 1 keV Ar+ 
ion beam, which was scanned over the target surface, at ^ = 0.0233 mA 
cm~^ with a fluence of 1.05 x lO'^^ ions cm~^ at three temperatures, namely 
-110° C, 23° C and 36° C. Under these conditions the pattern was only 
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Fig. 16. 500 X 500 nm^AFM images of a Si(lOO) substrate irradiated at normal 
incidence by 1.2 keV Ar"*" at different target temperatures: (a) 300 K, (b) 425 K. 
The vertical bars indicate 25 nm. PSD (c) and GID (d) curves for the patterns 
produced at 300 K (■) and 425 K (o). The vertical solid and dashed lines indicate 
the k = 2n/X (PSD) and g = 27r/A (GID) values for the 300 K and 425 K systems, 
respectively. 

obtained at 9 = 23°. In parallel, surface roughness increased markedly 
with temperature from 0.5 nm up to 76.8 nm. 

In addition, there is one study [60] where it is found that for low ion 
fluxes the dot size changes with increasing temperature according with the 
existence of Ehrlich-Schwocbel energy barriers [54, 150] whereas for high ion 
fluxes the dot size decreases with temperature. However, Si dots seem to 
arrange differently from other IBS nanodot structures induced on Si surfaces, 
speciflcally a clear pattern can not be visualized. 

Nanodot pattern dependence on ion flux 

For GaSb [57, 66] and InP [66] surfaces it was found that the pattern wave- 
length was independent of the ion current density or ion flux for the different 
ranges sampled, namely 10^^ - 4 x 10^^ cm~^ s~^ and 6.2 x 10^^ - 5 x 10^^ 
cm~^ s^^. respectively. For fixed InP targets irradiated at normal incidence 
by 1 keV Ar+ ions at fixed fluence and 23°, Tan and Wee [169] found that at 
low ion fluxes there was not any dot pattern but it appeared at jion = 0.0174 
mA cm~^. The pattern in-plane hexagonal order increased when ion flux was 
increased up to jion = 0.0233 mA cm^-^. In another study, where dot struc- 
tures did not form a clear pattern, it was proposed that for ion fluxes below 
220 /iA cm~^ Ehrlich-Schwoebel energy barriers [54, 150] dominate while for 
higher flux values the dot size decreased as ~ l/<?^/^ [60]. 
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In-plane order of nanodot pattern 

The nanodot patterns usually present short-range hexagonal in-plane order. 
The analysis usually employed for assessing their degree of order is made 
through the two-dimensional auto-correlation function C{r,t) of the AFM 
images 



which is a measure of how well a structure matches a space-shifted version of 
itself [193]. An example is shown in Fig. 10 where an AFM image of the dot 
pattern is displayed in panel (a) together with its corresponding 2-D auto- 
correlation (panel b). Here, six bright spots are clearly observed indicating 
the short-range hexagonal ordering of the dot pattern. 

When the symmetry of the pattern is large, it can be also assessed through 
the two-dimensional Fourier Transform (FFT) of the AFM images [197], al- 
though this is not the most common case, as IBS patterns usually give ringed 
FFTs [197]. 

Besides the symmetry of the pattern, a further issue is to quantify its 
degree of order in the pattern. As mentioned for the ripple patterns, two 
main approaches exist: (a) through size of the mean peak of the PSD of AFM 
images, and (b) using synchrotron techniques such as GISAXS and GID. 

Based on PSD data, Bobek and coworkers observed that the range of order 
of the nanopattern for GaSb surfaces increased appreciably after 100-200 s of 
Ar+ ion irradiation [12]. Also with this type of analysis the group of Frost 
found that the range of in-plane order decreased with increasing temperature, 
in the 268-335 K range, for rotating InP surfaces sputtered at ^ = 30° by 0.5 
kcV Ar+ ions [68]. 

Two studies have been performed for Si surfaces. For rotating Si substrates 
bombarded at ^ = 75° by either 0.5 keV Ar+ ions, by 1 keV Kr+ ions or by 1 
keV Xe+ ions [199]; Zibcri and coworkers obtained that the lateral correlation 
length, ^, increased with sputtering time up to a value close to 170 nm, in the 
window sampled, when Ar+ ions were employed. In contrast, for both Kr+ and 
Xe^, C saturated at ~ 145 nm, and ~ 120 nm, respectively. Saturation took 
place earlier for Xe+ . Moreover, the change of C, with ion energy depended on 
the ion species. Thus, for Kr+ and Xe+ the order increased with ion energy 
until the normalized correlation length, C/A, saturated for 1 keV. In contrast, 
for Ar"*" ions C, had a maximum value at 0.5 keV. 

The second system consists of a fixed Si target bombarded at normal inci- 
dence by 1.2 keV Ar+ ions. In this work, we used both AFM and synchrotron 
techniques to study the pattern dynamics in terms of coarsening and ordering 
as functions of sputtering time. In Fig. 17 we display the data obtained on 
Si by our group [72] and by the group of Frost [199] in terms of the normal- 
ized ordered domain size (i.e, the ratio of the lateral correlation length to 
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the pattern wavelength) versus ion dose. For the rotating Si target the order 
seems to increase with ion dose irrespective of the ion species, although satu- 
ration is observed for 1 keV Kr"*" ions but not for 0.5 keV Ar"'". In principle, 
it appears that the order of the pattern is larger for the rotating substrate 
configuration than for the fixed configuration when Ar"*" ions are employed. 
However, it should be noted that for 1 keV Ar"*" ions the order for the rotating 
configuration dropped to less than half that obtained for 0.5 keV Ar+ ions. 
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Fig. 17. Normalized ordered domains 
size, C/'^) "^s ion dose for: fixed Si(OOl) 
target sputtered by 1.2 keV Ar+ ions 
under normal incidence with ^ = 240 
/iA cm~^ (•); rotating Si target sput- 
tered by 0.5 keV Ar"*" ions (■) and 1 
keV Kr+ with $ = 300 ^A cm"^ (o). 
Data for the rotating substrates have 
been adapted from [199]. 



Regarding the quantification of the range of order in the pattern, SPM- 
based techniques pose a problem due to their locality. These surface character- 
ization techniques sample a relative small area of the surface and lack enough 
statistics to a certain extent. Thus, in a recent work [72], we have used AFM 
and synchrotron based techniques (namely, GISAXS and GID) to quantify the 
ordering of the same IBS patterns. The corresponding results for the normal- 
ized ordered domain size are presented in Fig. 18. Clearly, AFM data saturate 
at a value very close to 3 while GID data roughly increase monotonously in 
the sampled temporal window up to a value close to 10. This difference is 
due to the improved sampling statistics of GID with respect to AFM [72]. 
Thus, although AFM is a well suited technique for routine characterization 
of IBS nanopatterns, GID and GISAXS can provide us with more reliable 
quantitative data. 

Nanohole or nanopit patterning 

In principle, hole nanopatterns can also be induced by IBS. The fact that 
nanodots or nanoholes are produced on the target surface depends only on 
the anisotropies of the collision cascades [87] . After the ion sputtering theory 
by Sigmund [153] this shape depends only, for a given target material and ion 
species, on the ion energy. However, despite such a theoretical prediction, there 
is not clear evidence, up to now, of hole or pit pattern production by IBS. Two 
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Fig. 18. Normalized ordered domain 
size, (/\ obtained by AFM (•) and 
GID (o) versus sputtering time for a 
Si(OOl) target sputtered by 1.2 keV 
Ar"*" ions under normal incidence with 
# = 240 fxA cm"^ 



main pit-structures have been produced by IBS: the so-called cellular struc- 
tures [33, 37, 69, 84] and hole structures on semiconductor heterostructures 
[105]. The former look like network of hole-like structures, with a relatively 
wide distribution of hole sizes, resembling those obtained in plasma etching 
of silicon [192]. In contrast, those induced on semiconductor heterostructures 
display a more homogeneous hole size distribution with typical hole diameter 
of 170 ± 30 nm and hole spacing of 190 ± 40 nm. 

Recently, we have achieved the production of nano-hole patterns by IBS 
on silicon surfaces. These patterns display a characteristic wavelength similar 
to that obtained on the standard IBS experiments on dot pattern production. 
As an example wc present in Fig. 19a a typical AFM image of a nanohole IBS 
pattern. The hole structures are 2-3 nm deep and have a lateral size in the 
30-40 nm range. It should be noted that both the hole depth and lateral size 
can be unclerestimated because of tip convolution effects. The inset of Fig. 19a 
displays the auto-correlation of an AFM image taken over an area of 740 x 740 
nm^. Although there is some distortion, the short-range hexagonal order is 
clear. Finally, in panel (b) the radially averaged PSD function of image (a) is 
presented. Similarly to the results obtained above for nanodot patterns, the 
PSD presents a clear peak that is associated with the characteristic pattern 
wavelength that in this case is close to 50 nm. 

General considerations 

In general, the b{4iaviors found for ripple and dot nanopattern formation are 
analogous. Nevertheless, the ripple patterns present specific properties, due 
to the anisotropy of the process, such as ripple orientation, ripple propagation 
and shadowing. Analyzing the behaviors reported for the changes of both 
pattern morphologies with other parameters, the following map of trends can 
be drawn (we collect the main experimental observations in tables 3.1 and 
3.1): 

• IBS patterns are produced with different ion species, even with those im- 
plying reactive sputtering. 
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Fig. 19. (a) 1 X 1 pm^AFM image of a nanohole pattern induced on Si surfaces 
irradiated by 1 keV Ar"*" ions under normal incidence for a fluence of 4 x lO^'^ ions 

cm^^. Inset: two-dimensional auto-correlation function corresponding to a surface 
area of 740 x 740 nm^ in which the short-range hexagonal order is observed. The 
horizontal bar corresponds to 165 nm. (b) Radially averaged PSD function corre- 
sponding to the image of panel (a). The peak observed corresponds to an average 
inter-hole distance of fs 50 nm. 



• IBS patterns are produced on a large variety of substrates: amorphous, 
semiconductors and metals. 

• The pattern symmetry reflects that of the experimental set-up (i.e., 
anisotropic or isotropic). 

• In general, the pattern in-plane order increases with fluence (i.e, sputtering 
time). 

• Frequently, the pattern wavelength coarsens with time to finally saturate. 

• The surface roughness initially increases sharply to later on either saturate 
or increase at a lower rate. 

• At low temperatures, the pattern wavelength does not depend on target 
temperature and ion flux, but it increases with ion energy. 

• Ripple patterns: 

- Ripples run perpendicular or parallel to the ion beam projection direc- 
tion depending on the ion beam incidence angle. 

- Ripples propagate with a non-uniform velocity. 

- Shadowing effects appear in the ripple temporal evolution for long sput- 
tering times and large 6 values. 

- Ripples are produced by ions with energies ranging from just a few 
hundreds of eV up to 10^ eV. 

- At high temperatures, A increases with target temperature following 
the Arrhenius behavior and decreases with ion flux and ion energy. 

• Dot patterns: 

- Patterns are produced mainly with in-plane hexagonal order but also 
with square symmetry (only for oblique incidence and rotating targets). 

- To date, dot patterns are only produced by low-energy ions. 
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- Different dependences with target temperature have been found imply- 
ing cither no change, increase or decrease of A with increasing target 
temperature. In addition, for InP the pattern symmetry changes from 
hexagonal to square with increasing temperature. 

In general, the behaviors for ripple patterns seem to be better established 
than for dot patterns. This may be related to the fact that dot pattern produc- 
tion by IBS is relatively recent. In particular, the dependence of the morpho- 
logical properties of these patterns with temperature, ion energy and fluence 
has to be systematically addressed in order to obtain a more general picture 
of the nanopatterning process. On the other hand, both types of IBS patterns 
share many behaviors, which is consistent with the corresponding theoretical 
models. 

Also, there is a further experimental issue to be addressed: the possible 
influence of technical parameters on IBS pattern production. This important 
problem has been only studied by Ziberi and coworkers [194]. They obtained 
that the settings of the Kaufman ion-gun, the one that they employed, can af- 
fect the IBS nanopattern. In particular, they found that the divergence of the 
ion beam as well as the angular distribution of the ions within the ion beam 
influenced the pattern formation. If this happens already for a given ion-gun, 
we should likely expect some differences also when using other ion-gun types. 
Thus, it is convenient to perform systematic experiments on different targets 
and under different experimental conditions with the same equipment. This 
approach would allow a more direct comparison and contrast of the differ- 
ent experimental findings, which also would contribute to improve contrast 
between these results and the theoretical predictions. 



3.2 Pattern formation in single-crystal metals by ion beam 
sputtering 

The general observation of ripple and nanodot formation by off-normal and 
isotropic ion beam irradiation is only valid in the case of amorphous materials 
or for those that amorphize upon ion bombardment, such as single-crystal 
semiconductors. The case of metals should be treated through a different 
approach since the above properties do not occur. For example, ripples can be 
produced at normal ion incidence, or isotropic patterns may be formed under 
off-normal irradiation [146] . A large number of experimental results reported 
in metals has been compiled by Valbusa et al. [177]. 

The peculiar behavior of metal surfaces relies in the nature of the metallic 
bond. Due to its non-directional character, metals do not amorphize upon ion 
bombardment, at least for low fluence and up to energies of a few keV [177]. 
Therefore, the surface retains its properties, the ion bombardment being only 
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responsible for producing vacancies or vacancy aggregates at the surface [41], 
which increase the already large surface diffusion in these systems. Another 
relevant constraint for surface diffusion in metals comes from the presence of 
Ehrlich-Schwoebel (ES) energy barriers for adatoms to descend steps [54, 150]. 
These contributions add more complexity to surface nanostructuring by IBS 
in the case of metals. 

The diffusive regime [177] for pattern formation in metals appears when the 
ES energy barrier or a preferential diffusion path determines the final pattern 
characteristics. In this regime, the erosion process is limited to incorporating 
a larger number of diffusive particles (such as atoms or vacancies) that align 
along the more thermodynamically favorable directions. This means that, in 
contrast with amorphous or semiconductor materials, the intrinsic properties 
of the surface reflect onto the resulting surface pattern. For example, the 
pattern will reflect the intrinsic isotropic characteristics of the material as 
in the case reported for Pt(lll) by Michely et al. [115]. In the presence of 
diffusion anisotropy, as shown in Fig. 20 for Ag(lOO) [146], ripples may be 
observed even under normal ion incidence, or the ripple orientation may vary 
with temperature. In this context, temperature governs thermal activation of 
diffusion pathways, leading under certain conditions to isotropic or anisotropic 
configurations. Additional relevant conclusions extracted from Fig. 20 are that 
the pattern wavelength increases with temperature as a result of enhanced 
surface diffusion and that the pattern coarsens with sputtering time. It has 
also been observed that the pattern wavelength under this diffusive regime 
depends only slightly on the ion energy. 




Fig. 20. STM images (350 x 350 nm^) of Ag(llO) surfaces bombarded at normal 
incidence with 1 keV Ar+ for 15 min. at (a) 160 K, (b) 230 K, (c) 270 K, (d) 290 
K, (e) 320 K and (f) 350 K (taken from Ref. [177]). 
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The erosive regime in metals [177] is attained only for near glancing in- 
cidence angles and at low temperatures in order to inhibit thermal surface 
diffusion. In this regime the pattern formation can not be correlated with any 
symmetry of the surface, such as the orientation of the crystal. This fact is 
shown for Cu(llO) surfaces in Fig. 21 where the ion beam projection is mod- 
ified with respect to the surface orientation without relevant change on the 
pattern characteristics [147]. This implies that IBS patterns can be aligned 
in directions that are not thermodynamically favorable, this fact being one of 
the major potentials of surface nanostructuring by IBS in contrast with other 
techniques such as MBE. In this regime, the ripples are parallel to the beam 
direction, their wavelength depends linearly with the ion energy and both, 
wavelength and roughness, increases with fluence [177]. 




Fig. 21. STM images (400 x 400 nm^) of Cu(llO) surfaces bombarded at 70° with 

respect to the surface normal and 1 keV Ar^ for 15 min. at different angles between 
the beam projection and the surface orientation (adapted from Ref. [177]). 

Finally, ion bombardment of metals does not always results in the for- 
mation of a regular pattern on the surface. In this case, the surface may 
additionally display kinetic roughening or not depending on the conditions. 
An example of these results was obtained by Krim et al. [101] on Fe surfaces 
bombarded with 5 keV Ar+ at 6* = 25°. 

The different parameter windows for pattern formation, smoothening or 
roughening have been studied by Chason et al. [31] for Cu(OOl) surfaces as 
a function of ion dose and target temperature. In this scheme, the surface 
remains flat for low ion fluxes and high temperatures whereas it roughens 
for high fluxes and low temperatures. For intermediate temperatures (200- 
350 K) the diffusive regime dominates and, finally, high ion fluxes and high 
temperatures imply an erosive regime. 

3.3 Pattern formation in thin metal films by ion beam sputtering 

In the previous section, we have seen that ripples can be also produced by 

IBS on single-crystal metals for which ES energy barriers [54, 150] play such 
a crucial role, infiuencing the stability of the ripple morphology even at room 
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temperature [177]. However, a different scenario has to be considered for IBS 
nanostructuring of thin nictaUic films. These systems are usually polycrys- 
talline and mainly formed by grains that are randomly oriented, and with 
sizes that are usually much smaller than the film thickness [85, 180]. For such 
systems, the existence of ES barriers becomes improbable due to the lack of 
well-defined atomic steps at the surface [85] . These facts support the effective 
existence of an isotropic surface diffusivity rather than an anisotropic one. In 
this sense, thin metal films subject to IBS nanostructuring processes would 
be akin to the case of amorphous or amorphizable surfaces. Two groups have 
reported their findings for such systems. 

Karmakar and Ghosc have produced ripple structures on Co, Cu, Ag, Pt 
and Au thin film (thickness in the 30-200 nm range) surfaces irradiated at 80° 
by 16.7 keV Ar+ ions [90]. They found that the ripple wavelength as a function 
of the target element was qualitatively consistent with the behavior predicted 
by Makeev et al. [112]. Under this approach A depends mainly on the lateral 
and longitudinal straggling widths of the ion cascade, as well as on the angle 
of incidence. Karmakar and Ghose obtained these values from the computer 
code SRIM [201], and found qualitative agreement between the theoretical and 
experimental behaviors. However, there are sizeable quantitative differences. 
This is not new, quantitative disagreement between calculated and measured 
A values having been previously pointed out by different groups [93, 57, 73]. 
Another interesting finding refers to the stability of the induced ripple patterns 
at room temperature [90] in comparison with those produced on single crystals 
[177]. This fact is an indication that the thermally activated diffusion energy 
barriers in thin polycrystalline films are comparatively higher than in single 
crystal metal surfaces, which agrees with the assumption on the absence of 
ES barriers. In fact, measurements of activation energies for adatom surface 
diffusion on various polycrystalline materials are consistent with this result 
[142]. The studies of Stepanova and coworkers consisted in the irradiation 
by a 1.2 keV Ar+ ion beam of a 50 nm-thick Cu film deposited on glass 
or silicon substrates [162, 163]. When the irradiation angle was oblique, at 
= 82°, a clear ripple pattern was induced. These experiments allowed to 
discard any possible infiuence of the underlying substrate (either glass or 
silicon) in terms of composition and crystalline structure on the production 
of ripple patterns on the top metal thin film. Another set of experiments was 
performed under the same conditions but at normal incidence. In this 
self- assembled network of w 20 nm-sized ring-like Cu features was observed on 
a Si02 substrate after 1.2 keV Ar+ beam etching of the Cu/substrate interface 
[162]. 

The previous studies suggest that IBS of interfaces might be a reasonable 
process to fabricate metal nanopatterns on nonmetallic substrates. Also, the 
specific properties of these thin film polycrystalline metal systems make them 
suitable of study through to the same approach that is employed to understand 
IBS nanopatterning of amorphous/amorphizable surfaces. 
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4 Theoretical approaches 

As seen in the previous sections, typically the IBS induced nanopatterns fully 
evolve in macroscopic time and length scales (minutes and several microns, 
respectively). It is at these scales where, e.g., interaction among ripples can be 
seen to lead to order improvement with fluence or, rather, to eventual disorder 
in heights. Although detailed knowledge on the phenomenon of sputtering is 
rapidly and consistently developing through (microscopic) Moleciilar Dynamic 
type of studies (see e.g. [15, 131] and references therein, and Sec. 2), the scales 
that are currently reachable to these methods remain in the 1 ms and 50 nm 
ranges. Monte Carlo (MC) and continuum methods (CM) can probe larger 
scales so that these are the approaches that we will consider in what follows. 

4.1 Sigmund's theory of sputtering 

As mentioned in section 2.5, in a classic work Sigmund analyzed the kinetic 
transport theory of the sputtering process [153]. Assuming an infinite medium, 
he found that, in the elastic collision regime at the energies of a few keV 
where electronic stopping is not dominating, the deposited energy can be 
approximated by a Gaussian distribution near its maximum. Specifically, the 
density of energy spread out in the bulk by an ion with kinetic energy E is 
given by 



where the origin of the r' = {x' , y', z') coordinate system is placed at the 

impact point of the ion within the surface; z' is aligned along the ion beam 
direction, and x' and y' belong to the perpendicular plane to z', see Fig. (22). 




(5) 



In Eq. (5), = [(2^)3/^/^2] 



-1 



is a normalization constant and, due to the 




Fig. 22. Two-dimensional profile of the 

distribution of energy deposition, Eq. (5) . 
a and are the transversal and longitudi- 
nal widths of the Gaussian, respectively. 
The reference frames considered in the 
text are also shown. Adapted from [112]. 



initial momentum of the ion, the maximum of energy deposition occurs at 
a distance a along the ion trajectory inside the bulk. The longitudinal and 
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transversal straggling widths of the distribution are a and /i respectively. If 
mass difFcrcnces between substrate and ion are not too large, the quality of 
the approximation is reasonable and expression (5) provides a good approxi- 
mation for polycrystalline and amorphous targets [153]. Despite the fact that 
a depends on the microscopic details of the interactions between the ion and 
the bulk, at intermediate energies (10 — 100 keV), it is usually considered to 
be proportional to E. Within this range of energies, a, v and a are of the order 
of a few nanometers. This point will be relevant when information about the 
evolution of the topography is obtained. 

As we noted above, correlation of the sputtering yield with the surface 
geography is a crucial issue. In Ref. [154], Sigmund showed that the topogra- 
phy of the surface can indeed influence the magnitude of the rate of erosion 
and, provided an analytical description that describes the increase of yield for 
geometries different from the flat morphology. It is assumed that the speed of 
erosion at a point O on the surface is proportional to the amount of energy 
deposited there by the ions, with a proportionality constant A that is charac- 
teristic of the substrate and depends on the atomic density of the target n^, 
the atomic energy of connection in the surface \Jq, and a proportional con- 
stant Co which is related to the square of the effective radius of the potential 
of effective interaction according to 



Under these hypotheses, it is possible to obtain the mathematical expres- 
sion for the rate of volume eroded in O. This is simply given by 



The integral extends to the region where the impact of the ions contributes 
to energy deposition at O. The term <?(r') represents the local flux. We have 
specified its dependence with r' to show the corrections due to the local ge- 
ometry to the homogenous flux ^q. From (7), Sigmund derived the rate of 
erosion for diverse artificial geometries, such as a flat surface followed of an 
inclined plane or a vertex, in the case of normal incidence. 

The previous description, as already observed by Sigmund [154] , and briefly 
noted in Sec. 2.5 above, implies the occurrence of a morphological instability. 
Let us suppose that we irradiate a certain surface with an homogenous flow, 
as shown in Fig. 23. We can verify that distances OA and OB are smaller than 
O' A! and O'B' due to the geometry of the interface. This implies that, for this 
energy distribution, the penetrations of ions at A and B induce large energy 
deposition at O than the impacts on A! and B' at O' . As the rate of erosion is 
proportional to the deposited energy, erosion is faster at O than at O' . Thus, 
valleys are excavated more quickly that crests amplifying initial differences in 
heights. Sigmund suggested that an alternative process that flattens the sur- 
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face must exist and he proposed atomic migration as a mechanism to correct 
this instabiUty. 



The theory of Sigmund is the basis for most of the later continuum ap- 
proaches to the dynamics of surfaces undergoing by IBS. In principle, as a the- 
ory of sputtering, and as mentioned in Sec. 3.1, it is known to have limitations. 
Within the range of energies we are considering, two of the most conspicuous 
ones refer to: (i) the behavior of the sputtering yield for large incidence angles; 
(ii) crystalline targets. Related with the former, Sigmund's theory predicts an 
steady increase of the yield with incidence angle, not being able to account 
for the well-known decrease of S at glancing incidence due to more efficient 
projectile reflection and channeling effects at the surface. Nevertheless, recent 
continuum studies [35] have shown that improved distributions accounting for 
this effect (sec in [35] and [129]) do not modify qualitatively the morphological 
predictions to be derived below from the energy distribution (5). Similarly, 
channeling effects limit in principle the applicability of Sigmund's theory to 
crystalline solids [point (ii)] , by inducing sizeable departures from the energy 
distribution as described by (5). Thus, different distributions must be con- 
sidered, as done e.g. in [61] for the particular case of Cu targets. However, 
as seen in this work, again these modifications do not alter the qualitative 
morphological implications of the continuum theories to be described below, 
which justifies (at least for temperatures at which Ehrlich-Schwoebel barriers 
do not dominate surface diffusion [177]) the strong similarities in nanopattern 
formation by IBS on metals, as compared with amorphizable targets. 





|t>4| O'A' 
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As a morphological theory, again Sigmund's [154] has some limitations: 
it docs not predict the aheration of the morphology during the process for 
scales much larger than penetration depth for ions; as surface diffusion is 
not considered, the wavelength of the pattern needs to be of the order of 
the length scales appearing in the energy distribution Eq. (5); the effects of 
surface shadowing or redeposition are not considered; it does not predict the 
time evolution of the morphology and how it affects the rate of erosion. Thus, 
an additional physical mechanisms and a more detailed description of the 
surface height are needed in order to derive a morphological theory with an 
increased predictive power. 

4.2 Monte Carlo type models 

(Kinetic) MC methods have a long and successful tradition in the context of 
IBS [200, 201]. For our morphological purposes, they have found increased 
application during the last decade. Their main advantage is that, in princi- 
ple, they allow to make a more direct connection to microscopic relaxation 
mechanisms than e.g. continuum approaches. In addition, they do not require 
analytical approximations, lending themselves, rather, to numerical simula- 
tions that can reach reasonably large scales (around 10 s and around 100 
monolayers of erosion). However, as compared with CM, dependencies of ob- 
servables with physical parameters are less direct, accessible scales are shorter 
and universal (i.e. substrate/ion independent) properties are harder to assess 
unambiguously. Actually, the two types of approaches — MC and CM — to 
some extent complement each other so that, although our focus will be mostly 
on CM, we consider shortly some of the main results of the most recent MC 
studies. 

In order to access large scales as mentioned, MC methods consider physi- 
cal relaxation rules in which microscopic details — such as specific interactions 
among ions and target atoms — are coarse-grained out and replaced by effec- 
tive dynamical rules. These have to represent, on the one hand, the events 
leading to erosion by sputtering and, on the other hand, competing relaxation 
mechanisms such as surface diffusion. For the sputtering processes, starting 
with the study of [45] many MC models consider an erosion rule that is either a 
discretization [45] or, rather, a direct implementation [81, 173, 17, 187, 188, 32] 
of Sigmund's continuum law for energy deposition due to linear collision cas- 
cades (with corrections for large local incidence angles [161, 163]). Actually, 
this energy deposition rule has been extended even to MC studies of different 
physical processes, such as e.g. laser induced jet chemical etching [118]. Nev- 
ertheless, other MC models do consider more microscopic (albeit still approx- 
imate) rules for erosion, such as the binary collision approximation similar to 
that employed in SRIM simulations [200, 201], as done in the works by Koponen 
and collaborators [96, 97, 98, 99, 100]. Or else, still different approaches to the 
sputtering events consider dynamical creation of surface adatoms/advacancies 
[29, 126, 164]. 
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As for the additional surface relaxation mechanisms, most of the works 

mentioned restrict their nature to surface diffusion (there being additional 
physical possibilities, such as redeposition, viscous flow, etc.) for adatoms 
and/or advacancies. The rates for diffusion vary among the different studies, 
being either dependent on the energies of the initial and final states — within a 
local equilibrium approximation — or, rather, following temperature-activated 
Arrhenius laws; see a recent overview in [32]. 

Although there are some differences in the morphological dynamics pre- 
dicted by the various MC models, there is agreement in a number of ba- 
sic properties, such as: (i) the occurrence of a morphological instability at 
small/intermediate times, that is well described by CM; (ii) at oblique ion 
incidence, the non-linear evolution of the ripples as reflected in their non- 
uniform transverse motion; (iii) the saturation (stabilization) of the surface 
roughness at sufficiently long times. Differences among various models do oc- 
cur in their predictions: for instance, in some models the dominant pattern 
wavelength remains fixed in the course of time [45, 96, 97, 98, 99, 100, 32], 
while in other models a consistent wavelength coarsening phenomenon oc- 
curs [81, 162, 187, 188, 31, 163] with properties (like the coarsening exponent 
value) that depend on physical parameters like the substrate temperature. 
Then, there are properties such as the the kinetic roughening of the inter- 
face height fluctuations at long times — about which some of these studies do 
not extract conclusions (if they do not focus on such large scale issues), while 
other works reach different quantitative results (regarding e.g. the values of 
the scaling exponents, although most predict their occurrence). 

4.3 Continuum descriptions 

The wide separation in scales between the microscopic events that take place 
in the target during the irradiation process, and the large-scale response of its 
surface morphology recall a similar between microscopic and collective motion 
in the case of e.g. fluid dynamics. Thus, it is natural to expect [46] that some 
continuum description in the spirit of, say, the Navier-Stokes equations might 
be appropriate in the case of nanopattern formation by IBS. The advantage of 
continuum descriptions, if available, is that they provide compact descriptions 
of complex physical phenomena. Moreover, these frameworks are frequently 
more efficient computationally for the study of large scale properties, whose 
generic properties can be faithfully described by them [102]. Nevertheless, 
for the type of non-equilibrium phenomena which are our present focus, the 
relevant continuum equations are most times non-linear, their study requiring 
either analytical approximations or numerical simulations. 

Dynamics of the surface height 

In order to put forward a continuum desc;ription of these non-equilibrium sys- 
tems, Sigmund's theory [153, 154] constitutes an ideal starting point. This is 
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so because Eq. (7) already provides an explicit expression for the most rel- 
evant kinetic quantity describing the sputtering problem, namely, the local 
velocity of erosion. In turn, one may add different contributions to the right 
hand side of this equation if additional physical mechanisms contribute to the 
local variation of the target height, and this has indeed been done in order to 
improve the physical description. Still, for theoretical work, (7) is technically 
a very complex expression. The insight of Bradley and Harper (BH) [13] is 
that, under reasonable physical approximations, (7) becomes a closed time 
{t) evolution equation for a single physical field, the target height h{x,y,t) 
above point {x, y) on a reference plane. It is within such (simple) single-field 
approach that most of the recent continuum descriptions have circumscribed 
themselves. However, we will see in this section that this program encounters 
consistency problems — that question its physical applicability — unless the 
dynamics of additional independent fields arc taken into account. The evolu- 
tion of these fields — e.g. the density of species that diffuse on the surface — 
also modifies the local height velocity. This does not mean that a description 
by a closed equation for the height field is inappropriate but, rather, that to 
correct such description, one has to necessarily take into account for its deriva- 
tion physical mechanisms that complement Sigmund's formula. In this section 
wc recall the main results obtained along BH's pioneering "single-field" height 
approach, while the next section will (mostly) consider recent developments 
that incorporate the dynamics of additional physical fields to the continuum 
description. 

Bradley and Harper's theory 

Bradley and Harper (BH) [13] derived a partial differential equation to de- 
scribe the evolution of the morphology assuming that the variation in the 
height of the surface is smooth when seen at a scale that is comparable with 
the average penetration length of the ions. By expanding the target height 
as a function of the surface geometry up to linear order in the surface curva- 
tures, they obtained a closed evolution equation for the surface height. In the 
following we recall the main steps of the derivation. 

Define a local coordinate system with the z axis oriented along the surface 
normal at O and the x and y axes lying on the perpendicular plane to z (i.e. 
in the tangent plane through O) where x is aligned along the projection of 
the ion beam onto this plane as shown in Fig. 22. 

Assuming that the function which represents the height of the surface, 
z{x,y)^ is single valued and varies smoothly, we can approximate the height 
in the neighborhoods of O as 



where i?,,. and Ry are the principal radii of curvature. The simple formula (8) 
is enough to our purposes for two different reasons. First, for smooth surfaces 
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the average penetration length a is smaller than the radii of curvature and 

only ions penetrating the target a distance smaller than a from the point 
O contribute appreciably to erosion there (since energy deposition decays 
with distance as a Gaussian). Second, it is assumed that the surface principal 
directions — along which curvature is maximized or minimized — are aligned 
along the x and y axes. This holds at least for the cases in which structures 
are aligned either parallel or perpendicular to the projection of the incident 
beam (which is the usual case, as observed in previous sections). 

We can obtain the local speed of erosion at O using Eq. (5). To this end, 
one needs to change coordinates in Eq. (5) to the {x,y) frame in Fig. 22 
and take into account that the local flux of ions at any point of the surface, 
r = {x,y,z), is related to the homogenous flux, (^o, by #(r) = #ocos7(r), 
where 7(r) is the angle between the beam direction and the local normal 
direction. By Taylor expanding to linear order in the curvature radii and 
using Eq. (7), BH obtained an expression for the erosion rate at O that reads 



Vo = So 



(9) 



where the functions Fq, and Fy only depend on (I>q, E, A, the local angle of 
incidence 70 and the distances, a, a and fi of Eq. (5). 

Once the dependence of the erosion rate with the local morphology has 
been obtained, we can obtain a continuum equation for the local surface 
height, Z = h{X, Y), in the laboratory system of reference (X, Y, Z) sketched 
in Fig. 22. Thus, we define the component Z to be perpendicular to the initial 
surface, X as parallel to the projection of the ion beam onto the initial flat 
surface, and as Y perpendicular to X and Z. Within our linear approximation 
(see [112]) we can assume dth ~ —Vq and = 70 + dxh, where 9 is the angle 
between the ion beam and the normal to the uneroded surface, and using (9) 
one obtains [13] 

dh dh d^h d^h 

g^ = -VO + 7.^ + + r^ygy,, (10) 

with 



vo=SoFo{0), j^ = SodeFo{e), = SoF,{e), Uy = SoFy{e), (11) 

where any nonlinear terms have been neglected and we have used = 
—d\h and = —dyh. Here, So and Fq are always positive and, therefore, 
the surface height decreases at constant speed SqFq. 

In order to get information about the evolution of the surface, we will 
assume the following perturbation with spatial frequencies kx and ky, and 
amplitude A at time t: 



h{X,Y,t) = ^e'('=^^+'=^^)+'"*. 



(12) 
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If we substitute this expression into (10), we obtain the real part of the dis- 
persion relation, w. This reads 

7?.ea;(k) = —Vxk\ — Vyky, (13) 

and the imaginary part is 

Imw(k) = 7a;fcx- (14) 

The imaginary part of w(k) indicates the magnitude and direction of the 
speed of transverse in-plane motion of the disturbances. Since Iq is an in- 
creasing function oi 9 [13], the perturbation moves in the —X direction (that 
is, upstream with respect to the projection of the beam). On the other hand, 
TZe,u){\s.) describes the rate at which the amplitude of a perturbation with 
wave- vector k grows or decays with time. Given that at small angles of inci- 
dence [13] Vx{d) < and VyiO) < 0, any spatial perturbations on the initial 
surface grow exponentially in time. Since V^/i is positive at the bottoms of 
the valleys, at these points the rate of erosion is larger than at the top of 
the crests where V^/i < 0. Therefore, the negative signs of and Vy are 
the mathematical expression of Sigmund's morphological instability. At these 
small incidence angles, moreover, iyx{(^) < Vy{^) and perturbations grow faster 
along the X axis than along the Y. Hence, at these angles, the ripples crests are 
perpendicular to the projection of the ion beam onto the surface as observed 
experimentally. For the critical angle 9c, one precisely has Vxi^c) = i^y{(^c) 
while for larger angles VyiO) < VxiO) and the ripple crests align with the ion 
beam projection. 

These results predict some of most basic experimental featiires of IBS 
ripples but, without an additional mechanism which stabilizes the system, 
disturbances of arbitrarily small length-scales would increase exponentially 
without bound. On the other hand, from Eq. (10) one would expect the ripple 
wavelengths to be of the order of the distances involved in the description, 
that is to say, of the order of the penetration length a. This does not happen, 
rather, the experimental ripple wavelengths are frequently almost two orders of 
magnitude larger than a. In order to solve these problems, BH incorporate the 
effects of the surface diffusion of thermal origin by introducing an analogous 
term to that derived by MuUins [121] for isotropic surfaces. With this aim, a 
term —BV^V^h is included into (10) to obtain the following Hnear equation 
which describes the evolution of h 



dh dh d^h d^h 

VQ + Ix^TTT + t^x^TTTK + Vy 



dt~ " ' '"ax ' "ax2 ' *^^ay2 



- BV^V^/i. (15) 



Proceeding as in Eqs. (12)-(14), the real part of the dispersion relation 
modifies into 

7^ew(k) = -Uxk% - Uyk^ - Bk^, (16) 

that is negative for wave vectors of sufficiently large k (i.e. short length scales). 
However, there is a band of unstable long-wavelength modes for which (16) 
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is positive. The orientation and wavelength of the observed ripple structure 
corresponds to that mode which maximizes (16); this is 27r(2B/i/^)i/2 if the 
ripples are oriented along the X axis as a result of Vx < or 2iT{2B/vyY/'^ 
if the ripples are oriented along the Y axis as a consequence of Uy < v^- 
With the aim of obtaining an order of magnitude estimation of this quantity, 
Bradley and Harper deduced that for normal incidence, a sputtering yield of 
2 atoms per ion, a ~ cr ~ /x ~ 10 nm, and a typical value B ~ 2 x 10~^^ cm^ 
s^^, the value of the linear wavelength is A; ~ 5 ^m, in reasonable agreement 
with the typical distance between ripples observed experimentally. 

If surface diffusion is thermally activated, as effectively occurs for high 
temperatures and small flows, it is easy to obtain how the wavelength of the 
pattern varies with the temperature, T. Since v^^y do not depend on T and as- 
suming that B verifies B ~ (1/T) exp(— Z\_E/fcBT), where AE is the energetic 
barrier to activate surface diffusion and fc^ is the Boltzman's constant, the lin- 
ear wavelength A; must verify A; ~ (1/T^/^) exp(— Z\£'/2A;Br). On the other 
hand, it is also possible to obtain a relationship between ripple wavelength 
and the kinetic energy and flux of ions, E and <?0: respectively. Assuming that 
a, a and v are proportional to E and independent of one gets A; ~ E~^l'^ 

— 1/2 

and A; ~ ' . 

Working still within a linear approximation to the surface height in a sim- 
ilar spirit to the introduction of surface diffusion in (15), some alternative 
physical mechanisms have been proposed, specially in order to account for 
the lack of pattern (ripple) formation at either low temperatures [28] or small 
angles of incidence [19]. Thus, a term of the form —F\k\h{k,i) added to the 
(Fourier transformed) right hand side of Eq. (15) can describe the effect of 
viscous flow in the bulk onto the surface dynamics, with F being a coefficient 
that depends on the surface free energy and bulk viscosity [121, 122]. For low 
enough temperatures at which thermal surface diffusion is hampered, bulk 
viscous flow could dominate the surface dynamics even to the extent of pre- 
venting ripple formation for oblique angle incidence [28] or dot formation for 
rotating substrates [40] . Alternatively, the fact that a fraction of the sputtered 
atoms that move close to the surface are recaptured has led to arguing [19] 
for the effective description of such a process by inclusion of a term with the 
form +fc'^eccC?^/i, where I'recc would be a positive, angle dependent coefficient. 
We will see later that also direct knock-on sputtering events would have effec- 
tively an analogous continuum description, thus corresponding to stabilizing 
mechanisms that can counterbalance Sigmund's instability. 

Higher order corrections 

As we have seen, BH's linear theory already predicts many important fea- 
tures of the nanopatterning process, but leaves out a number of additional 
properties. For instance, it cannot predict saturation of the ripple amplitude 
or their non-uniform lateral motion. In order to account for these features. 
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one needs to generalize the BH theory, the most natural procedure being the 
extension of their perturbativc expansion in surface derivatives in (5)- (7) up 
to higher order. In this process one needs to approximate the surface height 
locally through a Taylor expansion, such as BH did, only that the order of 
approximation is higher. Along the way, not only non-linear contributions do 
arise, but also corrections appear to previous linear terms. These thus provide 
contributions of erosive origin to linear relaxation mechanisms such as trans- 
verse pattern motion, and surface diffusion. This program has been carried 
out in [44, 111, 112]. The result is the following equation (the coefficients A 
with various indices appearing in this and the following height equations are 
not to be confused with the pattern wavelength) 

dh__^^ , ^'^^ I j7 dhd^h dhd'^h 

dt dx ^ dx^ ^ dxy'^ ^ dx dx'^ ^ dx dy"^ 

d'^h d^h A^ /dhV fdhV , , 

"""OxW dx^ ^"9^4' 

where all coefficients depend on the experimental parameters from Sigmund's 

theory [112], like in BH's equation (10), which now becomes a linear low order 
approximation of Eq. (17). In the first line of (17), and except for the constant 
average velocity vq, all terms provide linear and non-linear contributions to 
the transverse ripple motion. The second line reflects the dependence of the 
sputtering yield with the local curvatures (as in BH's equation) and slopes, 
while the third line contains (linear) surface diffusion terms that include both 
thermal (as in BH) and erosive contributions. This type of contributions break 
in general the x — > —x symmetry while respecting symmetry under y — > —y, 
dependence of all parameters on the incidence angle being such that for normal 
incidence rotational in-plane symmetry is restored, Eq. (17) becoming the 
celebrated Kuramoto-Sivashinsky (KS) equation: 

^^=-vo + i^V'h - DV'h + ^{Vhf, (18) 

where we have used that, for ^ = 0, D^cx = Dyy = Dxy/2, = i^y = 
As = Ay = Ao, = = iy = ^\ = ^1 = ^ [112]- 

The main morphological implications of Eqs. (17) and (18) have been dis- 
cussed elsewhere [112]. Here, we briefly mention the main conclusions. First 
of all, from the point of view of these equations, the behavior predicted by 
BH appears as a short-time transient, and all BH predictions — on the de- 
pendence of e.g. the ripple structure with incidence angle — still carry over 
here. Moreover, for temperatures at which surface diffusion is activated, the 
dependences of the ripple wavelength with energy, temperature and flux are 
as in BH, at least for short to intermediate times. If T is low enough that the 
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only contribution to surface diffusion is from erosive origin, these dependen- 
cies arc modified and compare better with experiments at low temperatures, 
see [112] and Sec. 3 above. For intermediate to long times, the nonlinear terms 
in Eqs. (17), (18) are such that the exponential growth of the ripple ampli- 
tude is stabilized, yielding to a much slower power-law increase of e.g. the 
surface roughness with time. At normal incidence, the surface displays kinetic 
roughening for long time and length scales, and this is also the case at oblique 
incidence, as long as the nonlinearitics \x.y{d) have the same signs. Other- 
wise, Eq. (17) features cancellation modes as first identified in the anisotropic 
KS equation, that is the following particular case of Eq. (17) in which the 
propagative terms with coefficients S^x.y and i7i 2 (being arguably irrelevant 
to the asymptotic limit) are simply neglected [44, 144]: 




Cancellation modes are height Fourier modes with wave vector in the unstable 

band, for which the non-linear terms cancel each other, leaving the system 
non-linearly unstable, and inducing ripples which are oriented in an oblique 
direction that is parallel neither to the x nor to the y axis [144]. 

Eq. (18) has, on the other hand, the capability (beyond linear equations 
such as that of BH) of predicting both "dot" and "hole" production, depend- 
ing on the sign of the non- linear parameter A, which in turn depends on 
the characteristics of Sigmund's distribution [87]. However, conspicuous non- 
linear features still remain beyond description by the KS equation (18) and its 
generalizations (17), (19). First, the "dot" or ripple structures described by 
these non-linear equations are characterized by a dominant wavelength that 
remains fixed in time and does not coarsen for any parameter values, so that 
experiments in which coarsening occurs cannot be accounted for. A stronger 
limitation is that the patterns described by (17), (18), (19) disorder in heights 
to the extent that the size of ordered domains of ripples or dots essentially 
restricts to a single structure. There is thus no lateral ordering of the dots or 
ripples. Actually, the KS equation is known as a paradigm of spatio-temporal 
chaos in the wider field of Non-Linear Science [119]. 

Following the program sketched at the beginning of this paragraph, a nat- 
ural step is to carry on still further the perturbative study of Sigmund's local 
velocity of erosion pioneered by BH. At the present stage, one could close Eq. 
(18) (we now consider normal incidence for simplicity) by including non-linear 
terms that are quadratic in the height field and fourth order in space deriva- 
tives, reaching an equation of the form (in reference frame comoving with the 
eroded surface) 



(20) 
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Indeed, such a higher order generahzation has been performed [94]. However, 
the expressions of Ai and A2 in terms of Sigmimd's parameters arc such that 
the coefRcients of the two nonhnearities in Eq. (20) have the same signs for 
all physical parameter values. Unfortunately, as shown shortly after [26, 95], 
this introduces cancellation modes that seriously question the mathematical 
validity of (20) for our physical system: writing the equation in Fourier space, 
we get 



where FT denotes Fourier Transform. For the unstable mode ko (that indeed 
occurs physically, see the experiments of IBS of Pd(OOl) in [94, 26]) such that 
^0 = A1/A2, the nonlinear terms cancel each other, so that the amplitude 
of this mode explodes exponentially. We are seemingly left with a matter-of- 
principle limitation, namely, the simpler theoretical approach that is limited 
to studying Sigmund's local velocity of erosion (at sufBciently high linear and 
non-linear orders) meets mathematical limitations before being able to cover 
for the various experimental features of the patterns we wish to study. 

In the following paragraph we will take a wider viewpoint in which the 
dynamics is more complete, in the sense that the surface height will be cou- 
pled to an additional physical field describing the flux of adsorbed material 
that diffuses on the near-surface layer. This procedure will be seen to provide 
an improved description solving some of the above physical and mathematical 
shortcomings. Before doing so, we note that additional studies exist that focus 
on the evolution of the height field only. Thus, the standard result in the field 
of Pattern Formation (see references in [136]) that a linear damping in the 
KS equation induces the appearance of ordered patterns, has directly led to 
the proposition in [59] of (a modified) damped KS equation to describe the 
experiments described in Sec. 3. This modified equation actually becomes the 
standard damped KS equation after an appropriate non-local time transfor- 
mation [178]. Although the natural anisotropic generalization has been duly 
proposed [179], these equations unfortunately do not allow much improvement 
in the continuum description of nanopatterning by IBS for several reasons: (i) 
there is no connection to phenomenological parameters (the equations are not 
derived from any model but are, rather, argued for on a phenomenological 
basis); (ii) damped generalizations of the KS equation are known (see [27] 
and references therein) not to allow for wavelength coarsening for any param- 
eter values, which leaves out of these descriptions many of the experimentally 
observed patterns; (iii) similarly, at long times the fiuctuations of the PSD 
function predicted by any of the damped generalizations of the KS equation 
are cut-off and do not show the type of power-law like behavior that is seen, 
even in experiments in which there is no coarsening [195, 199]. In contrast, a 
more successful extension of the BH-type approach has been its generalization 
in [35] to surfaces with steep slopes. Recall that starting with BH's, all con- 
tinuum approaches mentioned (in contrast e.g. with MC models) work within 



9/t(k) 
dt 



{-pk^ - Dk^)h{k) + (Ai - X2k^)FT [{Vhf] , 



(21) 
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a small-slope approximation that allows neglection of higher order terms in 
Eqs. such as (17)-(20). Remarkably, a suitable generalization for arbitrary 
slope values has been seen in [35] to lead to a non-linear equation whose trav- 
eling wave (shock) solutions compare well with experiments on the motion of 
the walls of pits excavated by a FIB. These results possibly provide important 
clues for a more complete model of nanopatterning by IBS. 

Coupling to diffusive surface species 

Due to physical and mathematical considerations, the previous section leaves 
us with the need to enlarge the continuum description of nanopatterning by 
IBS. The expectation is that, by incorporating the dynamics of additional 
physically relevant fields, the effective height equation to be eventually de- 
rived improves its formal properties and its predictive power. Perhaps we 
could compare the situation with related fields such as e.g. the growth of thin 
films by physical or chemical vapor deposition (CVD) techniques. In CVD, 
for instance (see references in [46]), the standard continuum description arises 
precisely from the coupling between the local (growth) velocity and the dy- 
namics of the concentration field of the diffusing species that eventually will 
stick to the growing aggregate. 

There have been various attempts in the IBS context to combine the sur- 
face dynamics as predicted by a BH-type equation, with the evolution of 
relevant surface species, such as adatoms, addimers and surface traps (see 
[53, 30] and references therein). However, in this approach no explicit feed- 
back mechanism is provided from the dynamics of such species onto the local 
variation of the surface height so that the dynamics cannot be described by a 
closed system of equations. On the other extreme, there is a recent proposal 
in which a full Navier-Stokes (thus, highly coupled) formulation is proposed 
to describe ripple transverse motion onto a glass surface [2]. Trying to reach 
a balance between complexity and completeness in the physical description, 
one can seek for a formulation that, while simpler than a full hydrodynamic 
model, still provides an explicit coupling between the surface topography and 
the evolution of the relevant diffusive fields. For the case of ripple dynam- 
ics in the different (macroscopic!) context of aeolian sand dunes [42], such is 
the spirit of the so-called "hydrodynamic" approach, in which one sets up a 
system of coupled equations that describe the height of the eroded substrate 
profile, h, and the thickness of a mobile surface layer, R. Although there are 
relevant differences between both physical systems — the size of the structures 
is roughly seven orders of magnitude larger in aeolian ripple formation than in 
IBS systems — , the similarity between the mechanisms of diffusion and erosion 
and the analogous form and behavior of the patterns in both systems seem to 
suggest that these processes could be modeled by similar formalisms. 

This program has been followed in Refs.[6, 7, 25, 123], where the next 
equations are proposed to describe the evolution of the fields h and R: 
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dtR = (1 - -r'"^-V -3, 

dth = -r^^ + r'"^. 



(22) 
(23) 



In (22) and (23), n''(R,h) and r'"^{R,h) are, respectively, rates of atom 
excavation from and addition to the immobile bulk, (1 — (/)) = (j) measures 
the fraction of eroded atoms that become mobile, and the third term in Eq. 
(22) describes motion of mobile material onto the surface. In this way, local 
redeposition is allowed ii (f) ^ [104], while the viscous near-surface layer R 
is provided with a dynamics of its own [175]. 

In Refs. [6, 7], a linear dependence of F^^ and F"''^ with the local geometry 
of the surface is considered to study the linear stability of this system. These 
studies reveal that depending on the values of the parameters, the ripple ori- 
entation could be aligned in any target direction, one of the limitations of this 
model being that the model is not well defined in the absence of redeposition 
^ = 0. 

More detailed mechanisms of erosion and addition are explicited in Refs. 

[25, 123], where nonlinear terms for the rate of excavation arc introduced. 
These additional mechanisms are seen below to induce richer pattern dynamics 
than previous models for the surface height only. Unlike Aste and Valbusa's 
and acolian sand ripples models, in [123] the momentum in the direction of 
the projection of the beam directly transmitted from ions to superficial atoms 
is neglected. Rather, a diffusive term for mass transport onto the surface is 
considered that, in the case of isotropic amorphous materials, is given by 
J = —DVR, where D is a thermally activated constant. Another feature 
of [25, 123] is to consider a non-zero amount of mobile material, Req, even 
in the absence of excavation (/^«^ = 0) o redeposition {(J) = 0) which could 
be thermally induced. This term allows us to write the rate of addition in 
similar form to the Gibbs-Thompson expression for surface relaxation via 
evaporation-condensation, namely. 



where 70 is the mean nucleation rate for a flat surface (on the xy plane) and 

72x) 72y ^ describe variations of the nucleation rate with surface curvatures 
(and are positive if we assume that nucleation events are more likely in surface 
valleys than in protrusions). Note that, in (24), the full thickness of the mobile 
layer is affected by the shape of the surface. 

In [123], a expression similar to the non-linear generalization (17) of BH's 
local erosion velocity is considered for the rate at which material is sputtered 
from the bulk. If the beam direction is in the xz plane, we have, 



-Tex = Q!o [1 + ctixdxh + a2xdlh + a2ydyh + as^id^h)'^ + asy{dyKf 

-{dxh){a4xdlh + a4ydyh)], (25) 



where parameters reflect the dependence of F^x on the local shape of the 
surface and are function of the angle of incidence, ion and substrate species. 



(24) 
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ion fliix, energy, and other experimental parameters just as the coefficients of 



Oblique incidence 

In a weakly nonlinear analysis of this system, the length-scale separation we 
have already mentioned between the time-scales associated with erosive and 
diffusive events is seen to play a cnicial role. Thus, in particular the diffusive 
field R relaxes much faster than the target height h, making it possible within 
a multiple-scale formulation to solve perturbatively for the dynamics of the 
former and derive a closed effective equation for the evolution of the latter 
which, in the case of oblique incidence [123, 125], reads 



where these new coefficients are related to those of Eqs. (22)-(25) [123]. 

As in the nonlinear continuous theories shown in the previous section, 
the reflection symmetry in x is broken in Eq. (26), but it is preserved in 
the y axis. Equation (26) generalizes the height equations in [13, 112, 137] 
within BH approach to IBS; the main difference between (17) and (26) is 

(2) 

that additional terms A)^- appear in the later. These terms result important 
to correctly describe the evolution of the pattern as we will see bellow. To 
the best of our knowledge, Eq. (26) is new, and indeed has a rich parameter 
space. Numerical integration within linear regime retrieves all features of the 
ripple structure as predicted by the BH theory, i.e., dependence of the ripple 
wavelength with linear terms, and ripple orientation as a function of 9. 

Top views of the morphologies described by Eq. (26) at different times are 
shown in Fig. 24 for different values of the parameters. In all these examples, 
the ripples arc increasing their size in the course of time while disordering in 
heights for long distances, whereas the form of the ripples can vary appreciably 
depending on parameter values. In some cases, as in the first and third rows 
of Fig. 24, the ripples are disordered longitudinally, whereas in the second 
row, the ripples are longer and more straight. In the third row, the ripples are 
also very disordered in height and tend to group themselves in domains which 
contain, approximately, three ripples each. The ripple coarsening seen in Fig. 

(2) 

24 actually requires the presence of Xl J , whose magnitude and mathematically 
correct sign are due to describing redeposition by means of the additional field 
R. When the values of these coefficients increase relative to A^-"'^'' , coarsening 
stops later, and the amplitude and wavelength of the pattern also increase, 
analogous to the result for one-dimensional interfaces studied in [124]. The 
coarsening exponent n will take an effective value that will be larger the later 
coarsening stops, and may depend on (simulation) parameter values [123]. 



Eq. (17). 




i=x,y 



+ J2 [-IC^jdfd]h + xfj^dKdjhr],{26) 
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(a) (b) (c) 




(d) (e) (f) 




Fig. 24. Top views of morphologies described by (26) at t = 10 (first column), 
t = 106 (second column) and t = 953 (third column) with 7^, = —0.1, fx = 1, 
= 0.1, = 1, ny = 0.5, = 0.1, Ai'^ = 1, A^'^ = 5, A,'2j = -5, K-ij = 1 (a), 
(b) and (c). (d), (e) and (f) with the same parameters as in (a), (b) and (c) except 
Ai'^ = 0.1. (g), (h) and (i) with 7, = 0.1, = 1, i^y = -0.95, f2, = -0.5, = 0.1, 
Ai^^ = 0.1, Xi^^ = 1.0, X'l^l = -0.5, A^^^ = -5.0 y ICi,j = 1. The x axis is oriented 
vertically and the y axis horizontally. 



Asymmetry and transverse in-plane motion of the pattern depends on 
terms with an odd number of derivatives in Eq. (26). Whereas the linear 
terms tend to make the ripples move with a constant velocity, the nonlinear 
term S,i accounts for non-uniform movement of the ripples as observed in some 
experiments [79, 2]. In Fig. 25 we show the evolution of a transversal section 
of the surface as described by Eq. (26), where the non- uniform movement and 
the asymmetry of ripples can be appreciated 
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Fig. 25. Cross sections at j/ = L/2 
of the two-dimensional surface given by 
Eq. (26) for times between t = Q and 
t = 1500 separated by a time interval 

A^ = 20 with = -2, Qi ^ 0, u^, = 1, 



vy = 0.1, ii = 3.5, A, 



(1) 



1, A. 



(1) 
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An example of a comparison of Eq. (26) with experiments is shown in 
figure 26 where coarsening of the ripples is clearly appreciated. As we see, Eq. 
(26) indeed captures essential properties of the evolution of the experimental 
topography. 




Fig. 26. (a), (b) 1x1 /im^ AFM views of fused silica targets irradiated under 
oblique incidence after 10 and 60 min. respectively; (c), (d) morphologies obtained 
by numerical integration of Eq. (26) at times t\ and 6ti. Panels (a) and (b) are 
taken from [63] with permission. 
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Normal incidence 

In the case that the surface is bombarded perpendicularly, the in-plane asym- 
metry introduced by the oblique beam disappears and, for materials that do 
not show crystallographically privileged directions, the height equation which 
describes the evolution of h reads [25, 123] 

dth = -uV^h - ICV^h + A(i) {Vhf + A(2) v2(V/i)2, (27) 

where, again, the coefficients arc related to the Sigmund's parameters. A 
stochastic generalization of this equation has been also proposed in the context 
of amorphous thin film growth [140, 141]. 

An important property of (27) is that, in the absence of rcdcpositcd ma- 
terial {(j) = 0), A^^' and A*^^-* have the same signs, precisely as occurred in the 
case of the purely erosive approach [94] discussed in the previous section. This 
makes Eq. (27) nonlincarly unstable [141, 26]. Hence, the description with a 
single height field seems to feature an intrinsic problem which is solved, if a 
fraction of redeposited material is large enough, through the two-field "hydro- 
dynamic" approach. 

In order to reduce the number of parameters and simplify the analysis of 
(27), we rescale x and y, t, and h by {K/vY/'^,K/v'^ and i^/Ai, respectively, 
resulting into a single-parameter equation as done in the one-dimensional 
counterpart of (27) studied in [124]. This reads 

dth = -V^/i - V^/i + {Vhf - rV'^{Vhf, (28) 

where r = —{v\'^'^^)/{KX^^'>) is a positive parameter which allows us to per- 
form a numerical analysis of the complete parameter space of Eq. (27) through 
an appropriate rescaling. 

In Figs. 27 and 28 the evolution of the height profile and corresponding 
auto-correlation functions are shown for r = 5. Starting from an initial ran- 
dom distribution, a periodic surface structure with a wavelength of about the 
maximum of the linear dispersion relation arises and the amplitude of h in- 
creases. Then, a coarsening process occurs and dots grow in width and height, 
the total number of them decreasing. The apparent coarsening is quantified in 
the plot of \{t) shown in Fig. 29 where the saturation of ripple wavelength at 
long times can be observed. Simultaneously, with dot coarsening, the pattern 
increase its in-plane order and eventually leading to an hexagonal dot array 
as shown in Fig. 28. 

As an example of comparison of this equation with specific experiments, in 
this paragraph we show the results recently obtained in [72] where an array of 
nanodots is obtained over a silicon substrate. In this work, the time evolution 
of nanodots is considered on Si(OOl) and Si(lll) targets irradiated at normal 
incidence with Ar+ ions of 1.2 keV. In figures 30 are shown characteristic 
AFM images of Si(OOl) and Si(lll) surfaces eroded at different times. The 
pattern is similar in both experiments: in both cases dots of 5-7 nm height 
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Fig. 29. Temporal evolution of the 
pattern wavelength. A, given by Eq. 
(28) for different values of r. For r = 
50, the dashed line represents the fit 
to a power law A ^ t"'**^. 



and 40-60 nm width are formed and group into short-range hexagonal order. 
For times larger than 20 minutes, the surface disorders in heights for larger 
distances (of the order of 6 nm in height and about 500 nm of lateral distance) . 
This long-range disorder (kinetic roughening) increases in the course of the 
time. In this experiment, the rate of erosion was determined experimentally 
in all the processed samples by partially masking them during the sputtering 
process and measuring the resulting step edge height with a profilomctcr. It 
was observed that the sputtering rate (SR) was 10% higher for Si(lll) than 
for Si(OOl). In order to understand the differences found experimentally be- 
tween the two surface orientations, we have integrated Eq. (27) numerically 
considering this fact. Thus, as the model neglects target crystallinity, we have 
used this relative difference in SR to simulate the IBS pattern evolution in two 
model systems. Insets of Figs. 30(a) and 30(b) display the two-dimensional 2D 
top- view images from the simulations for surfaces with a given SR at two dif- 
ferent simulation times (time and length units are arbitrary) . Similarly, insets 
of Figs. 30(c) and 30(d) correspond to the same simulations for a surface with 
a 10% higher SR, thus representing the Si(lll) orientation. The simulations 
reproduce in both cases the experimental coarsening and occurrence of a long- 
wavelength corrugation. In order to quantify the pattern features, the PSD 
functions corresponding to simulations were systematically analyzed. In Fig. 
9 is shown that the pattern wavelength saturates earlier for the surface with 
higher SR, but the lower SR case attains a larger final dot size. In addition, it 
is also shown that the correlation length also saturates earlier for a higher SR 
surface. Remarkably, both AFM and simulation results agree in estimating 
ordered domains to contain roughly three nanodots. This agreement between 
simulations and experimental observations allows us to conclude e.g. that the 
experimental differences observed between the pattern evolution on Si(lll) 
and Si(OOl) surfaces are due to their different SRs. 
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Fig. 30. AFM 3x3 fim^ images of 
Si(OOl) and Si(lll) surfaces sput- 
tered with 1.2 keV Ar"*" at normal 
incidence for different times (see 
labels). Insets: 2D views (length 
scales in arbitrary units) from nu- 
merical integration of Eq. (27) for 
different simulation times and sur- 
faces with low {u — 0.68, K. — 1, 
X^^^ = 0.028, and A'^' = 0.136) 
and high {ly = 0.75, /C = 1, A^^' = 
0.031, and A^^^ = 0.15) erosion 
rates. Reprinted from [72] with per- 



mission. 



Rotating substrate 

Following the ideas of Ref. [14], in Ref. [123] we obtained the evolution equa- 
tion for the surface in case the substrate is rotating simultaneously while 
bombarding. This reads 



dth = -i/,v2/i~/CrV*/i+A(i)(V/i)2-hAp)v2(V/i)2+A(=') V-[{V'^h) V/i] , (29) 



where a relationship between these parameters and those in Eq. (26) exists. 

Preliminar simulations of Eq. (29) [125] show that this equation also 
presents interrupted coarsening and an ordered array of dots, in which the 

(3) 

effect in the morphology evolution of the new term Ar is similar to that of 



Comparison between continuum models 

The continuum description of nanopattern formation by IBS is currently still 
open and making progress. As seen in this Section, to date there are several 
alternative descriptions available that share some of their predictions, while 
differing in several other aspects. It is also true that theoretical work is not 
even complete yet — in terms of analyzing systematically the dependences of 
morphological properties with physical parameters such as e.g. temperature 
or ion energy — due to two main reasons: (i) the main interfacial equations 
are nonlinear, thus not easily amenable to analytical solutions, while some of 
the physically interesting features such as the stationary pattern wavelength 
depend crucially on non-linear effects; (ii) the parameter spaces of these mod- 
els are frequently large, particularly in the oblique-incidence case. Actually, 
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some of the interface equations that we have been considering seem to be new 
in the wider contexts of Non-equiHbrium Systems and Non-Uncar Dynamics 
[124, 27, 46] even to the extent of providing examples of thus far unknown 
behaviors for problems of high current interest, such as coarsening phenomena 
[139]. 

As a partial summary, and for the sake of comparison, we have collected 
in table 4.3 the main available morphological predictions of the models that 
we have considered in a larger detail in the text. After comparison with the 
analogous experimental tables 3.1 and 3.1, the two-field model leading to Eqs. 
(26), (27) and (29) seems to date the theoretical description that can account 
for a larger range of experimental behaviors from within a single framework. 
Naturally, more work is still needed in order to reach a still more complete 
description of these physical systems with a larger quantitative predictive 
power. 



5 Applications of IBS patterned surfaces 

As mentioned in the Introduction, surface nanopatterning is the subject of 
intense research driven by the current road to miniaturization and the broad 
range of eventual applications. Therefore, the formation of self-organized pat- 
terns on surfaces by IBS is relevant from a fundamental point of view but also 
from a practical point of view. Whereas the understanding of the patterning 
process has been studied deeply from the theoretical and experimental points 
of view, the applications of these patters have been scarcely addressed and a 
further technological effort is required for their implementation in real devices. 
Here, we summarize some of the applications launched in the literature. 

One of the strengths of the IBS nanopatterning process comes from its 
universality since, as shown in section 3, it is applicable to metallic, semicon- 
ducting or insulating surfaces. However, in comparison with other techniques, 
the major advantages of this technique rely on its simplicity, high output 
(fast) and scalable patterned area. Among the pattern characteristics, we can 
highlight the tunable wavelength and presence of ordering. 

The first applications of IBS patterns were based on the production of 
micro-scale ripple patterns. In this case, the usage of these patterns for optical 
interfca-encc^ gratings was addressed [86]. Recently, optical activation of the 
surface due to the pattern formation has also been reported [128] which, apart 
from the surface morphology, may be also linked directly to the ion induced 
defects. 

Perhaps the most important and direct application of IBS patterning relies 
with the method based on the deposition of a given material layer on top of the 
target material intended to be patterned. This layer is then ion sputtered, this 
process inducing a nanopattern on top of it that, by continuing the sputtering 
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Table 3. Summary of main morphological predictions for some of the continuum models discussed in the text. NM stands for Not 
measurable because there is not a well-defined saturation value of the roughness. 
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process, is eventually transferred onto the target substrate. This strategy was 
previously used to generate a Ni dot pattern with a typical wavelength of 
250-300 nm and dot height of 13-15 nm [151]. In this case, the top-layer 
was a copolymer with self-organized polymer domains. These domains did 
show different sputtering rates when ion sputtered. In this case, these Ni 
nanopatterns were used to separate DNA molecules. 

The interest on IBS patterns was greatly increased by the report of nan- 
odot pattern formation by Facsko ct al. [55] as the distances involved were 
clearly in the nanometer range. Therefore, the potential application for quan- 
tum dot (QD) arrays fabrication [56] triggered further research. This route 
for QD fabrication is shown schematically in Fig. 31 and was experimentally 
proved by the successful photoemission of confined nanostructures [55]. De- 
spite this result, there are issues to be improved for the final applicability, such 
as the crystallinity of the nanostructures. since ion bombardment induces con- 
siderable amorphization on the surface. This constraint may not play a big 
role if the nanodot volume is mostly crystalline. In any case, re-crystallization 
or damage reduction may be attained by post-annealing treatments. 



1) Ion irradiation 2) Pattern formation 3) Pattern transfer 




Fig. 31. Sketch of the fabrication of quantum dot arrays by transferring the self- 
organized nanodot pattern induced by IBS to the interface of a semiconductor het- 
erostructure (adapted from Ref [56]). 

The IBS patterns can be also used as molds or masters to transfer or repli- 
cate the pattern to other surfaces with high functionality. This several-step 
process may be required when the surface can not be processed directly, such 
as the case of polymer materials, or if the desired pattern can not be obtained 
in a simple way (for example, when the diffusive regime is operating in met- 
als). Also, this imprinting path may be interesting for reducing production 
costs by faster processing or avoiding the use of vacuum environments. The 
applicability of IBS patterned surfaces for molding and replication has been 
shown by Azzaroni et al. [8] in polymer materials but the method has been 
successfully extended to metals [10] and even hard ceramic materials [5]. The 
replication method is shown schematically in Fig. 32a and consists of several 
sequential steps: (1) production of the IBS dot pattern on a silicon surface; 
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(2) formation of a self-assembled monolayer of octadecyltrichlorosilane (OTS) 
on top of the nanostructured surface; (3) subsequent growth of the material 
to be patterned. This step has to be performed at mild conditions (i.e, low 
temperatures or under non-high energetic sputtering conditions) in order to 
preserve the OTS layer; (4) thanks to the antisticking properties of the OTS 
layer the film deposited on top of it can be easily mechanically detached. Also, 
due to the extremely low thickness of the OTS layer the film surface in con- 
tact with it follows the surface morphology of the underlying nanostructured 
surface. Thus, the surface morphology of the detached film is the negative of 
the IBS nanostructured silicon surface. In Fig. 32b we show the molded sur- 
face of a diamond-like carbon (DLC) film grown by radio frequency sputtering 
following this method. The negative hole structure of the IBS dot surface is 
observed while the two-dimensional auto-correlation function (inset) shows 
that the short-range hexagonal order is still preserved. 




Fig. 32. (a) Replication method for IBS nanodot pattern transfer consisting in 4 
steps: (1) production of the IBS pattern on Si surfaces; (2) functionalization of the 
nanostructured Si surface by a OTS monolayer; (3) film deposition; (4) mechani- 
cal film detachment, (b) 750 x 750 nm'^ AFM image of a DLC film surface molded 
following the method described in (a). Inset: Two-dimensional auto-correlation cor- 
responding to an area of 545 x 545 nm^. The horizontal bar corresponds to 63 nm. 

Apart from the direct exploitation of the properties induced in nanopat- 
terned surfaces, much effort has also been devoted to the use of these surfaces 
as templates for further processing. For example, ripple patterns have been 
tested as substrates for alignment and manipulation of carbon nanotubes [74] . 
In this sense, one possible application would be to extend FIB patterning as 
a template technique for patterned growth of carbon nanotubes [36] to IBS 
nanopatterning. Also, the large range of surface roughness and characteristic 
wavelengths that are attainable by IBS patterns open a broad field of bioap- 
plications. For example, surface adhesion of biological entities (proteins, DNA, 
. . .) may be controlled for sensing, diagnostic or biocompatiability enhance- 
ment. Moreover, the nanopit patterns could find applications as supports for 
catalysis [148]. Another potential application of IBS patterns is for the pro- 
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duction of magnetic nanostructures [172]. Moroni and coworkers proved that 
the formation of nanoscalc ripples by IBS on Co films deeply affect their mag- 
netic properties [120]. Also, IBS patterns have been used as growth templates 
for the production of magnetic nanostructures. For example, it is well known 
that the magnetic properties of a superlatticc formed by magnetic materi- 
als depend on the final thickness of the stack. Therefore, irregularities in the 
thickness due to non conformal growth, i.e. preferential growth on the cavi- 
ties (ripples or holes) or on the top of the hills, onto a patterned surface can 
be used to define magnetic domains [34]. Another proposal is the use of the 
surface topography to induce shadowed deposition under oblique precursors 
incidence [171]. These routes are displayed schematically in Fig. 33. 




Fig. 33. Proposals for the formation of magnetic nanostructures using IBS pat- 
terned surfaces as growtli templates: (a) non-conformal growth (the different thick- 
ness result in patterned magnetic domains) and (b) shadowed deposition. 

Obviously, the use of patterned surfaces for many applications is not re- 
stricted to IBS patterns. However, the presence of the, mostly undesirable, 
amorphous layer is a peculiarity of IBS patterns that can be used for further 
applications. For example, the modulation of the amorphous layer in medium- 
energy induced ripples was very successfully exploited by Smirnov et al. [158] 
as selective doping implantation barrier. In this case, the production of rip- 
ples on silicon by nitrogen erosion leads to the formation of amorphous silicon 
nitride (SiN^) regions, as shown in Fig. 34. The difi^erent nitride layer thick- 
nesses can be used to obtain silicon regions with different degree of doping 
through, for instance, reactive ion-etching and As+ implantation processes. 

6 Open issues 

Despite the large amount of theoretical and experimental work described 
above, there are still issues that remain unclear and require further research 
efforts. From the experimental point of view, this is especially evident for the 
case of nanodot patterns due to their relative novelty. For the case of nanorip- 
ples, in spite of the abundant available data, there are also gaps to be filled. 
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Fig. 34. (a) Formation of a rijjple pattern and inherent modulation in the amor- 
phized layer as a result of ion bombardment at off-normal incidence, (b) Cross- 
sectional TEM image of a ripple pattern produced by N'^ ion bombardment of silicon 
and the formation of different SiN^; regions that are used as implantation barriers. 
Figure reprinted with permission from [158]. 



Perhaps, the least systematically investigated aspects concern transverse rip- 
ple propagation and the relevance of shadowing effects on the morphological 
evolution. The former is difficult to address due to its evident experimental 
complexity since in-situ, real time monitoring of the ripple morphology has to 
be carried out. The second issue, although quite specific, is interesting because 
it can provide important data for theory refinement, which requires exploring 
systematically the continuum models beyond the customary small-slope ap- 
proximation. With regard to nanodot patterns, there is still a certain lack of 
systematicity in the studies of the various targets that have proved amenable 
to such patterns when subject to IBS. For instance, the role of preferential 
sputtering on nanostructuring of hctero-scmiconductors has not been unam- 
biguously assessed yet. In addition, the effect of temperature on the nanodot 
morphology or the symmetry of their arrangement is not clear as different, 
and in some cases opposite, behaviors have been found. Another issue that 
will be interesting to address will be the systematic comparison of the nan- 
odot pattern properties for the same target material bombarded under normal 
incidence and under oblique incidence with simultaneous target rotation. In 
this sense, comparison is still incomplete between the IBS patterns produced 
with the same equipment on different materials. A related open issue is the 
influence of the characteristics of the ion gun optics on the pattern properties. 

Another line of research with important implications in the theoretical 
understanding of IBS patterning is that of pattern coarsening. Prom the data 
compiled, it is clear that many systems present coarsening but also there are 
others that do not. Thus, an important experimental investigation would be 
the study of the physical processes behind pattern coarsening. Once more, the 
results of these investigations would contribute to improving the theoretical 
models and hence deepen our understanding of IBS patterning formation. In 
this connection, an ambitious goal is to clarify the relationship between the 
more relevant parameters appearing in the continuum models with the specific 
properties of the target materials and ion species. This research would indeed 
help to understand quantitatively the different patterns properties found ex- 
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perimentally, such as for instance the different size and shape of dots found 
in Si, GaSb and InP. Moreover, it would be interesting to design specific ex- 
periments in order to verify alternative predictions from the various models. 
An example of this can be the experimental verification of recent proposals to 
produce novel nanopatterns under specific multi-beam irradiation conditions 
[22, 23, 24]. 

Finally, from a more applied point of view, the use of IBS nanostructur- 
ing processes for technological applications remains as an important challenge 
since, although several specific applications have been already developed, they 
have not yet proved clear advantages with respect to previously existing al- 
ternative routes. 
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